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by 
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SUMMARY 

To  assist  in  the  prediction  of  the  Earth-coverage  capabilities  of 
large  constellations  of  satellites,  studies  have  been  conducted  of 
patterns  based  on  regular  tesselations  of  a  plane  surface  and  of  others 
based  on  regular  polyhedra.  This  work  has  identified  trends  associated  with 
increasing  numbers  of  satellites  and  degrees  of  coverage,  permitting  extra¬ 
polation  from  the  known  capabilities  of  smaller  satellite  constellations. 

An  investigation  of  selection  criteria  has  emphasised  the  importance  of  the 
minimum  inter-satellite  distance,  which  provides  a  particularly  convenient 
criterion  for  use  in  preliminary  studies,  and  this  in  turn  has  been  used  in 
an  examination  of  some  non-uniform  orbital  patterns,  which  has  suggested  an 
improved  method  of  synthesising  such  patterns.  Results  obtained  have  been 
discussed  with  reference  to  the  requirements  of  the  GPS  navigation  satellite 
system. 
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I  INTRODUCTION 

Methods  of  achieving  continuous  coverage  of  the  whole  surface  of  the  Earth  by 
optimised  constellations  of  satellites  have  been  studied  by  various  organisations  for 
more  than  20  years,  this  in  support  of  a  variety  of  satellite  missions  such  as  communica¬ 
tions,  navigation  and  surveillance.  At  RAE,  the  author  has  made  a  number  of  studies  in 
this  field,  many  of  which  were  collected  in  two  Reports  issued  in  1970  and  1977.  This 
present  Report  sunmarises  further  studies  performed  during  the  period  1977-1982. 

The  1970  Report'  described  analyses  which  had  been  made,  using  hand  methods,  of  two 

types  of  circular  orbit  constellations  which  were  described  as  star  patterns  and  delta 

patterns  respectively;  by  the  time  it  was  issued,  work  had  started  on  a  computer  program 

2 

(COCO)  needed  for  further  analysis  of  the  more  complex  delta  patterns.  A  197]  paper 

summarised  chat  report,  including  corrected  results  in  a  few  cases  where  early  computer 

3  . 

runs  had  shown  their  need.  A  brief  paper  issued  in  1973  contained  summarised  results  of 
comprehensive  computer  studies  of  delta  patterns  containing  up  to  15  satellites  providing 
up  to  four-fold  continuous  coverage  of  the  whole  Earth's  surface,  and  described  the 
method  of  classifying  patterns  which  had  been  adopted  for  use  with  the  computer  program. 

4 

The  1977  Report  provided  fuller  details  of  these  studies,  and  of  later  extensions  to 
cover  delta  patterns  of  up  to  25  satellites  providing  up  to  seven-fold  continuous  coverage 
of  the  whole  Earth's  surface;  it  also  contained  an  analysis  of  the  Earth-tracks  followed 
by  delta  patterns,  including  a  method  of  identifying  those  patterns  which  would  follow  a 
single  repetitive  Earth-track  at  a  particular  orbital  period,  and  identified  certain 
series  of  patterns  which  both  followed  single  non-self-crossing  Earth-tracks  and  included 
many  of  the  patterns  providing  the  highest  standards  of  coverage.  This  Report  also 
described  the  program  COCO  used  for  coverage  studies. 

In  Refs  I  and  2  star  patterns  were  defined  as  consisting  of  multiple  orbits  (with 
equal  numbers  of  equally-spaced  satellites  in  each)  sharing  cwo  conmon  points  of  inter¬ 
section  (which  may  be  considered  as  the  poles  of  a  reference  plane),  with  all  ascending 
nodes  (with  respect  to  the  reference  plane)  lying  in  one  180°  arc  of  Che  refer¬ 
ence  plane  and  all  descending  nodes  in  the  other,  and  with  equal  intervals  between  adjac¬ 
ent  ascending  nodes;  delta  patterns  were  also  defined  as  consisting  of  multiple  orbits 
with  equal  numbers  of  equally-spaced  satellites  in  each,  but  with  all  orbits  having  the 
same  inclination  5  to  a  reference  plane  and  with  ascending  nodes  evenly  spaced  around  the 
reference  plane.  Later  reports  and  papers  concentrated  on  the  delta  patterns,  which  were 
found  to  give  generally  superior  results.  Ref  3  introduced  the  identification  code  T/P/F 
for  delta  paccerns  which  had  been  adopted  when  writing  the  computer  program  COCO;  here 
T  is  the  total  number  of  satellites  in  the  pattern,  P  is  the  number  of  equally-spaced 
orbital  planes  of  inclination  5  between  which  they  are  equally  divided,  and  F  is  a 
measure  of  the  phase  difference  between  satellites  in  adjacent  planes  such  that,  when  a 
satellite  in  one  plane  reaches  its  ascending  node,  one  of  the  satellites  in  the  adjacent 
plane  having  a  more  easterly  ascending  node  is  F  units  of  360°/T  past  its  ascending 
node.  For  a  delta  pattern  P  may  be  any  factor  of  T  ,  including  I  and  T  ,  and  F 
may  have  any  integer  value  from  0  to  (P  -  1). 


This  series  of  reports  and  papers  discussed  standards  of  coverage  provided  by  the 
various  satellite  patterns  in  terms  (using  nomenclature  adopted  in  Ref  4)  of  a  parameter 
Rn  ,  the  angle  subtended  at  the  centre  of  the  Earth  by  the  radius  of  the  circumcircle  of 
three  sub-satellite  points  which  encloses  (n  -  l)  other  sub-satellite  points.  It  was 
demonstrated  (and  an  alternative  form  of  demonstration  is  given  in  section  2.2  of  this 
Report)  that  the  centre  of  such  a  circumcircle  is  the  locally  least-favoured  point  for 
seeing  at  least  n  satellites  at  maximum  elevation  above  the  horizon,  system  require¬ 
ments  usually  being  written  in  terms  of  a  minimum  elevation  angle.  The  quantity  R 

max,  n 

(with  a  small  'm')  was  defined  as  the  instantaneous  value  of  the  largest  circumcircle  radius 

in  che  constellation  for  degree  of  coverage  n  ,  and  &  as  the  common  inclination  of  the 

orbital  planes  of  the  constellation  to  its  reference  plane.  The  internal  configuration 

of  the  constellation  would  vary  systematically  during  an  orbital  period,  and  \fax  n 

(with  a  capital  'M')  was  defined  as  the  largest  value  of  R  during  a  complete 

max ,  n 

orbital  period,  for  any  particular  inclination  5  .  R^^  n  should  preferably  be  as 

small  as  possible,  in  order  to  make  the  minimum  elevation  angle  as  large  as  possible 

(with  satellites  at  a  given  altitude)  or  to  permit  the  satellite  altitude  to  be  as  low  as 

possible  (for  a  given  minimum  elevation  angle);  by  varying  the  common  inclination  6  of 

the  orbital  planes  it  is  possible  to  find  an  inclination  6  at  which  R„  has  a 

opc  Max,n 

minimum  value,  defined  as  Rj^  n  (with  'MAX'  in  capitals). 

The  main  objective  was  generally  taken  to  be  the  identification  of  the  pattern 

which,  for  a  particular  value  of  T  (and  perhaps  also  of  P  ),  could  provide  the 

smallest  value  of  Rj^  n  •  However,  it  was  considered  that  the  minimum  value  of  the 

inter-satellite  distance  D  was  also  of  some  significance  (as  discussed  in  section  2.3), 

with  any  value  less  than  about  3°  probably  being  unacceptable,  and  zero  minimum 

separation  almost  certainly  so;  hence  values  of  D  were  also  calculated.  In  similar 

manner  to  R  ,  D  .  was  defined  as  the  smallest  instantaneous  value  of  the  inter- 

max^  min 

satellite  distance  occurring  in  the  constellation;  as  the  smallest  such  value 

during  a  full  orbital  period;  and  ■  with  corresponding  inclination  4  ^  D  ,  as  the 

largest  value  of  D^n  which  could  be  found  by  varying  the  inclination  5  .  In 

Refs  1  to  4  only  values  of  D  .  at  6  were  calculated  and  presented  with  che  values 

Min  opt 

of  n  <  later  in  1977,  however,  after  preparation  of  Ref  4,  values  of  D^j.  were 

investigated,  and  the  results  were  described  briefly  in  a  paper^  issued  in  1978  and  are 
presented  more  fully  here. 

An  appendix  to  Ref  4  discussed  the  relevance  of  its  contents  to  navigation  satellite 
systems,  such  as  the  US  GPS  system  (using  Havstar  satellites)  which  was  then  expected  Co 
be  a  three-plane  24-satellite  system.  In  1980  it  became  known  that  the  initial  GPS 
system  would  contain  only  18  active  satellites,  in  a  pattern  still  to  be  determined;  as 
1 8-satellite  delta  patterns  had  been  examined  in  1974-5  using  the  program  COCO,  but  only 
some  of  these  results  had  been  listed  in  Ref  4  and  conmented  on  in  its  Appendix  A,  che 
relevant  computer  printouts  were  re-examined  and  an  unpublished  RAE  note^  on  the  findings 
was  prepared.  Information  was  subsequently  received  on  US  studies  of  some  non-uniform 
18-satellice  patterns;  a  brief  further  study  of  such  non-uniform  patterns  was  made, 
resulting  in  another  unpublished  RAE  noce7. 
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The  present  wricer  had,  during  preparation  of  Ref  1,  considered  the  possibility  of 

obtaining  guidance  for  the  study  of  satellite  constellations  from  a  study  of  the  regular 

9 

polyhedra,  and  to  this  end  examined  Cricchlow’s  book  on  the  characteristics  of  regular 
and  semi-regular  polyhedra  and  tesselations .  At  that  time  no  information  directly 
relevant  to  the  study  of  satellite  constellations  was  identified,  though  it  was 
recognised  chat  a  regular  polyhedron  represents  the  most  'efficient'  possible  distribu¬ 
tion  of  the  relevant  number  of  vertices  (or  satellites),  which  cannot  be  matched  on  a 
worst-case  basis  by  a  real  satellite  constellation  because  of  its  internal  relative 
motion;  also  that  the  equal  edge-lengths  of  a  regular  polyhedron  represent  a  maximisation 
of  the  minimum  inter-vertex  distance,  which  might  point  to  an  aspect  worth  considering  in 
satellite  constellation  studies.  When  information  was  received  on  a  l’S  proposal,  by 

Q 

Draim  ,  thac  satellite  constellations  might  advantageously  be  designed  so  chat  the 
different  orbital  planes  would  be  parallel  to  the  faces  of  a  regular  polyhedron, 
with  the  satellite  phasing  also  based  on  the  characteristics  of  a  regular  polyhedron,  it 
was  felt  that  this  proposal  deserved  further  examination  in  association  with  a  further 
study  of  Critch low's  work. 

It  was  found  chat  most  of  Oraim's  suggested  constellations  required  more  chan  one 
inclination,  which  would  probably  be  unacceptable  in  practice;  however,  one  which  required 
only  one  inclination  (Draim's  octahedron  constellation,  with  five  satellites  in  each  of 
four  planes  of  54.736°  inclination)  was  found  to  be  identical  to  the  delta  pattern 
20/4/2  which  the  present  writer  had  examined  in  1975.  From  the  work  on  delta  patterns  it 
appeared  that  the  optimum  inclination  for  pattern  20/4/2  would  be  considerably  less  than 
54°,  and  chat  some  other  20-sacellice  delta  patterns  could  provide  better  coverage;  it 

was  therefore  concluded  that  Draim's  specific  proposal  did  not  lead  to  optimum  results. 

9 

However,  on  re-examination  of  Critchlow's  work  it  was  realised  that  the  characteristics 
of  regular  tesselations  (ie  regular  partitions  of  a  plane  surface)  might  give  more 
guidance  than  those  of  regular  polyhedra  (which  provide  regular  partitions  of  a  spherical 
surface),  and  a  study  of  this  aspect  made  during  1982  is  reported  here.  This  study  has 
(i)  identified  trends  with  respect  to  increases  in  the  number  of  satellites  (T)  and  the 
degree  of  coverage  (n)  which  appear  relevant  also  to  regular  polyhedra  and  to  practical 
satellite  constellations;  (ii)  provided  a  basis  for  normalising  values  of  n  and 

^MIN  so  as  Co  Pera^c  extrapolation  of  results  established  previously  to  higher  values 
of  T  and  n  ;  and  (iii)  provided  a  basis  for  an  examination  of  a  potential  combined 
criterion  which  demonstrated  the  relative  importance  of  jjq  *  suggesting  its  use  as 
sole  criterion  for  preliminary  studies  to  establish  a  short-list  of  suitable  patterns  to 
meet  any  particular  requirement. 

To  provide  a  more  logical  development  of  the  discussion  than  would  a  strictly 
chronological  description,  this  account  of  these  studies  begins  in  section  2  with  a  dis¬ 
cussion  of  chose  selection  criteria  which  the  writer  has  found  useful  in  general  coverage 
studies  (as  opposed  to  the  more  detailed  criteria  which  may  be  developed  for  analysis  of 
specific  satellite  system  requirements).  Section  3  describes  the  1982  analysis  of 
tesselations  and  regular  polyhedra,  and  section  4  includes  both  the  application  of  this 
analysis  co  practical  satellite  constellations  and  the  results  of  the  1977  studies  of 

v 
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.  Section  5  covers  the  1980  investigations  of  18-satellite  patterns  potentially 
suitable  for  GPS,  including  both  the  re-examination  of  1974-5  computer  results  and  the 
brief  study  of  characteristics  of  non-uniform  patterns.  Finally,  the  principal 
conclusions  are  summarised  in  section  8. 

2  REQUIREMENTS  AND  SELECTION  CRITERIA 

2. 1  General 

The  problem  with  which  we  are  concerned  is  chat  of  providing,  in  the  most  efficient 
manner,  an  array  of  satellites  which  can  between  them  provide  continuous  multiple 
(n-fold)  coverage  of  the  entire  surface  of  the  Earth.  Since,  to  provide  whole-Earth 
coverage,  the  satellites  must  necessarily  be  in  multiple  intersecting  orbits,  the  pattern 
which  the  satellites  form  relative  to  the  Earth  is  constantly  changing,  and  the  coverage 
requirements  must  be  met  by  all  the  configurations  which  the  satellite  constellation 
takes  up  -  in  other  words,  we  must  ensure  that  the  requirements  are  3till  met  under 
worst-case  conditions.  This  is,  therefore,  a  more  complex  problem  than  that  of  providing 
an  array  of  earch-stations  able  between  them  to  provide  continuous  multiple  sightings  of 
a  satellite  in  any  orbit  (above  a  certain  altitude)  about  the  Earth,  since  the  earth- 
stations  would  form  a  fixed  pattern.  Nevertheless,  we  can  obtain  some  useful  insights 
from  examination  of  chis  simpler  earth-station  problem,  and  of  an  even  simpler  problem 
regarding  plane  arrays  of  sensors,  as  we  shall  consider  later. 

2.2  Circumcircle  radius  and  effective  horizon 

Coverage  studies  have  generally  been  conducted  in  terms  of  the  radius  from  the  sub¬ 
satellite  point  to  the  effective  horizon,  this  radius  being  measured  by  the  angle  it 
subtends  at  the  centre  of  the  Earth.  The  'effective'  horizon  reflects  the  particular 
requirements,  depending  on  satellite  altitude  and  on  the  minimum  acceptable  elevation 
angle  from  a  point  on  the  Earth's  surface;  for  example,  for  satellites  in  12-hour  circular 
orbits  to  be  visible  at  a  minimum  elevation  angle  of  5°,  this  condition  would  be  met 
within  a  radius  of  71.2°  of  a  sub-satellite  point  (from  Fig  1  of  Ref  I). 

Fig  la  represents  the  instantaneous  positions  of  the  sub-satellite  points  of  three 
out  of  a  constellation  of  satellites  in  circular  orbits.  A  circle  whose  radius  represents 
the  effective  horizon  distance  has  been  drawn  (with  a  solid  line)  around  each  sub- 
satellite  point.  These  circles  divide  the  area  covered  by  this  figure  into  three  types 
of  element:  within  those  marked  I,  one  satellite  is  visible;  within  those  marked  2,  two 
satellites  are  visible;  while  within  that  marked  0,  no  satellite  is  visible.  If  the 
radius  of  the  circles  were  increased  to  allow  for  an  increased  horizon  range 
{eg  satellites  in  higher  orbits,  or  lower  minimum  elevation),  the  area  marked  0  would 
shrink,  as  shown  in  Fig  lb,  and  full  single  coverage  of  this  area  would  eventually  be 
achieved  when  the  three  circles  all  passed  through  the  point  +  at  its  centre;  this  is 
the  centre  of  the  circumcircle  (drawn  with  a  broken  line)  of  the  three  sub-satellite 
points,  and  the  radii  of  the  three  circles  would  then  be  equal  to  this  circumcircle 
radius.  Thus  the  centre  of  the  circumcircle  of  three  adjacent  sub-satellite  points 
represents,  for  that  part  of  the  overall  pattern,  the  worst  case  for  meeting  the  single 
coverage  requirements;  and  the  centre  of  the  largest  of  such  circumcircles  is  the  worst 
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case  for  che  whole  pattern.  Considering  the  variation  of  the  pattern  during  an  orbital 
period,  the  largest  value  of  the  largest  such  circumcircle  radius  (i;  j)  defines  the 

pattern's  single  coverage  capabilities;  it  is  desirable  that  this  value  should  be  as 
small  as  possible. 

Fig  1c  represents,  in  similar  fashion,  the  instantaneous  positions  of  four 
neighbouring  sub-satellite  points.  The  point  ♦  at  the  centre  of  the  area  within  which 
only  one  satellite  is  visible  is  the  centre  of  the  circumcircle  of  the  outer  three  sub¬ 
satellite  points,  which  encloses  che  fourth  sub-satellite  point;  if  the  radius  of  the 
horizon  circles  were  increased  to  equal  that  of  the  circumcircle,  full  double  coverage 
of  che  area  would  be  obcained.  In  general,  the  centre  of  the  circumcircle  of  any  three 
sub-satellite  points  is  che  locally  critical  point  for  degree  of  coverage  n  ,  where 
In  -  1)  is  cne  number  of  other  sub-satellite  points  enclosed  within  the  circumcircle. 

Fig  Id  illustrates  the  situation  when  che  circumcircle  of  any  three  sub-satellite 
points  passes  simultaneously  through  one  or  more  (in  this  case  two)  additional  sub¬ 
satellite  poincs.  In  Fig  Id,  as  in  Fig  lc,  one  other  sub-satellite  point  is  enclosed 
within  the  circumcircle,  so  that  the  centre  of  the  circumcircle  is  again  the  locally 
critical  point  for  double  coverage.  However,  under  these  circumstances  it  should  also 
be  considered  a  locally  critical  point  for  triple  and  quadruple  coverage,  since  minimal 
changes  in  the  positions  of  one  or  two  of  the  sub-satellite  points  lying  on  the  circum¬ 
circle  could  leave  two  or  three  sub-satellite  points,  instead  of  only  one,  enclosed 
within  the  circumcircle.  Hence  if  the  circumcircle  passes  through  x  sub-satellite 
points  (where  x  is  not  less  than  3),  and  encloses  y  other  sub-satellite  points,  its 
centre  is  the  locally  critical  point  for  degrees  of  coverage  (n)  from  (y  +  I)  to 
(x  +  y  -  2)  -  though  since  x  will  often  be  equal  to  3,  these  values  will  often  be  equal. 

The  smaller  the  circumcircle  radius,  che  higher  is  the  minimum  elevation  angle  at 
which  the  required  number  of  satellites  (at  a  given  orbital  altitude)  can  be  seen  from 
any  point  within  the  circumcircle.  Likewise,  the  smaller  che  circumcircle  radius  che 
lower  is  the  orbital  altitude  necessary  for  che  required  number  of  satellites  to  be 
visible  (above  a  given  minimum  elevation  angle)  from  any  point  within  the  circumcircle. 

Thus  in  general  terms  it  is  desirable  that  che  largest  value  of  circumcircle  radius  (R^ax  n 
or  should  be  as  small  as  possible.  However,  once  the  requirements  for  a  particular 

system  have  been  chosen,  and  che  satellite  altitude  and  minimum  elevation  angle  (deter¬ 
mining  the  maximum  circumcircle  radius),  che  required  degree  of  coverage  (n)  and  the 
number  of  satellites  in  the  system  have  all  been  specified,  it  might  be  argued  that  it  is 
only  necessary  to  confirm  that  all  satellite  patterns  placed  on  the  short-list  for 
selection  meet  the  stated  requirement  for  maximum  circumcircle  radius,  and  that  the 
choice  between  them  should  depend  not  on  which  of  them  has  the  smallest  such  radius  but 
on  some  other  criterion. 

2.3  Inter-satellite  distance 

Another  criterion  which  has  previously  been  used  by  the  author,  with  increasing 
emphasis  in  successive  papers  (though  still  as  a  criterion  additional  to  the  circumcircle 
radius),  is  the  minimum  inter-satellite  distance.  In  the  earliest  paper'  it  was  simply 
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seated  chat  "The  minimum  satellite  separation  ...  may  be  of  interest  if  pairs  of 

satellites  are  to  be  used  for  position-fixing  ...”,  and  a  similar  comment  was  made  in 

Ref  2.  In  Ref  3  lists  were  provided  of  patterns  giving  minimum  values  of  R^^^  n  f°r 

different  values  of  T  (from  5  to  15)  and  of  n  (from  1  to  4) ,  with  the  inclination 

ac  which  chis  occurred  and  the  value  of  D...  at  this  inclination,  with  the  comment  that 

Min 

"Apart  from  considerations  of  interference,  [  the  value  of  Dj^n  ]  is  likely  to  be 

particularly  significant  where  position  determination  using  more  than  one  satellite  is 

required;  no  pattern  giving  a  value  of  Dw.  less  than  3°  has  been  included  in  the 

win 

table,  and  where  the  value  of  Dw.  is  small  an  alternative  pattern  giving  a  larger 

Min 

value  i3  listed".  Ref  4  said  "It  is  assumed  that,  as  a  secondary  objective,  the  minimum 
separation  between  any  two  satellites  in  the  pattern  should  be  as  large  as  possible.  The 
direcc  importance  of  this  objective  may  vary  according  to  the  system  application;  in  a 
satellite  navigation  system,  accuracy  may  well  increase  as  the  minimum  distance  between 
the  satellites  providing  the  fix  increases,  while  in  a  satellite  communications  system  it 
may  only  be  necessary  that  the  minimum  distance  should  exceed  some  fixed  value,  to  ensure 
that  interference  between  transmissions  in  the  same  frequency  band  is  acceptably  small. 
However,  it  also  has  some  indirect  importance  in  relation  to  the  main  objective;  the 
larger  the  minimum  distance  between  satellites,  the  more  uniform  the  distribution  of 
satellites  over  the  Earth's  surface,  and  hence  the  more  likely  that  the  pattern  will 
provide  relatively  favourable  values  of  the  maximum  distance  to  the  nth  nearest  sub¬ 
satellite  point,  for  all  relevant  levels  of  coverage".  Ref  4  later  noted  that  "For  most 
purposes  it  would  be  unsatisfactory  to  choose  a  satellite  pattern  in  which  pairs  of 
satellites  passed  very  close  to  one  another;  chis  ought  cause  radio  interference,  reduce 
the  number  of  independent  observations  available  from  the  system,  or  cause  other 
undesirable  effects,  even  if  the  possibility  of  physical  collision  were  discounted"; 
and,  referring  to  tables  presenting,  as  in  Ref  3,  values  of  n  with  associated 

values  of  4  ^  and  ^Min  *  *  commented  that  "If  it  were  important,  for  a  particular 

system,  that  the  minimum  satellite  separation  should  be  as  large  as  possible,  then  it 
would  be  appropriate  to  optimise  for  rather  than  for  Rj^  n  »  or  at  least  to  aim 

for  a  compromise  between  the  two". 

Ref  5  first  presented  some  results  of  optimising  for  »  and  showed  that  the 

largest  values  of  DMIN  were  associated  with  particular  sets  of  delta  patterns  which  had 
been  identified  in  Ref  4  as  being  associated  with  single  non-self-crossing  repetitive 
Earth-tracks .  It  noted  that  "two  types  of  requirement  are  likely  to  be  generally 
applicable  to  systems  requiring  continuous  multiple  whole-Earth  coverage-.  One  concerns 
the  level  of  coverage  required  (single,  double  ...  n-fold)  and  the  corresponding  minimum 
elevation  angle  to  the  nth  nearest  satellite,  and  the  other  the  minimum  angular  separa¬ 
tion  between  adjacent  satellites  in  the  pattern,  as  measured  at  the  centre  of  the 
Earth".  Hence,  over  the  series  of  papers,  the  inter-satellite  distance  has  been  treated 
initially  as  a  feature  which  might  be  of  interest  in  some  circumstances;  later  as  a 
secondary  objective  for  optimisation;  and  most  recently  as  the  object  of  one  out  of  a 
pair  of  requirements.  Hence,  if  the  circumcircle  radius  were  to  be  treated  only  as 
being  subject  to  an  upper  limit,  rather  than  as  the  subject  of  an  optimisation,  the 
inter-satellite  distance  might  be  a  suitable  choice  as  an  alternative  subject  for 


optimisation.  However,  there  are  circumstances  in  which  the  requirement  tor  Che  incer- 
satellite  distance  itself  might  only  be  that  it  should  exceed  some  specified  minimum 
value. 

2 .  -  The  possibility  of  a  combined  criterion 

Ref  4  suggested  the  possibility  chat  one  might  "aim  for  a  compromise  between  the 
two".  Certainly  there  are  difficulties  in  having  two  separate,  and  sometimes  conflict¬ 
ing,  criteria  for  optimisation;  later  in  this  Report  we  shall  examine  further  a  possible 
method  of  combining  the  two  objectives  into  one,  and  the  consequences  of  so  doing. 

However,  while  it  has  appeared  clear  that  patterns  having  a  large  value  of  the  minimum 
inter-satellite  distance  would  also  provide  favourable  values  of  circumcircle  radius, 
the  precise  nature  of  the  relationship  between  these  criteria  for  any  particular  pattern 
was  not  obvious;  and  it  was  thought  that  an  examination  of  the  simpler  earth-station 
problem,  using  regular  polyhedral  patterns,  might  help  to  clarify  this  situation. 

3  ANALYSIS  OF  REGULAR  AND  SEMI-REGULAR  POLYHEDRA  AMD  TESSELATIONS 

3 .  !  General  considerations 

It  seems  clear  that  an  ideal  solution  to  the  earth-station  problem,  that  of  most 
economically  providing  an  array  of  stations  able  to  crack  at  all  times  a  satellite  in  any 
orbit  exceeding  a  given  altitude,  would  be  provided  by  placing  stations  at  positions 
corresponding  to  the  vertices  of  one  of  the  regular  polyhedra,  if  the  altitude  were  such 
as  to  make  such  a  number  of  stations  appropriate.  However,  as  is  well  known,  there  are 
only  five  regular  (Platonic)  polyhedra,  namely  the  tetrahedron,  octahedron,  cube, 
icosahedron  and  dodecahedron,  having  respectively  4,  6,  8,  12  and  20  vertices,  so  the 
chances  of  meeting  any  specific  requirements  by  means  of  such  an  ideal  solution  are 
limited. 

9 

CritchLow  has  provided  a  comprehensive  review  of  the  regular  and  semi-regular 
polyhedra  and  tesselations ,  showing  that,  in  addition  to  che  five  regular  polyhedra, 
there  are  13  semi-regular  (Archimedean)  polyhedra;  of  these,  two  have  12  vertices,  four 
have  24,  one  has  30,  one  has  48,  four  have  60  and  one  has  120.  Thus  these  add  relatively 
little  to  the  chances  of  obtaining  an  ideal  solution  to  the  earth-station  problem.  They 
share  with  the  regular  polyhedra  che  property  of  having  only  a  single  value  of  edge-length 
lor  inter-station  distance),  but  whereas  the  regular  polyhedra  each  have  only  a  single 
face  configuration  consisting  of  a  regular  polygon  (triangular  for  tetrahedron, 
octahedron  and  icosahedron,  square  for  the  cube  and  pentagonal  for  the  dodecahedron),  the 
semi-regular  polyhedra  each  have  either  two  or  three  different  types  of  regular  polygons 
forming  their  faces. 

However,  we  may  appropriately  consider  the  regular  and  semi-regular  cesselacions  of 
a  plane  as  forming  extensions  of  these  series  of  regular  and  semi-regular  polyhedra, 
since  che  plane  may  represent  part  of  the  surface  of  a  sphere  of  infinite  radius. 

Critchlow  shows  that  (using  his  nomenclature)  there  are  only  three  regular  tesselations 
‘triangular,  square  and  hexagonal),  each  having  a  single  edge-length,  vertex  configuration 
and  regular  polygonal  face  shape;  eight  'semi-regular'  tesselations,  each  having  a  single 
edge-lengch  and  vertex  configuration  but  two  or  three  different  regular  polygonal  face 
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shapes;  and  14  ' demi-regular '  tesselations,  each  having  a  single  edge-length  buc  two  or 
three  different  vertex  configurations  and  from  cwo  to  four  different  regular  polygonal 
face  shapes.  Wich  such  tesselations  one  may  choose  to  consider  an  area  containing  any 
arbitrary  number  of  vertices,  from  zero  to  infinity. 

Whereas  we  have  previously  discussed  the  optimisation  problem  in  cerns  of 
'satellites'  or 'ear th-s cations ' ,  these  are  inappropriate  terms  to  use  in  reiati-n  to  a 
plane  surface,  so  we  shall  speak  instead  of  arrays  of  'sensors'  deployed  at  positions 
corresponding  to  the  vertices  of  tesselations.  For  example,  we  may  imagine  an  array  of 
sensors  deployed  on  a  flat  plain  to  track  the  movements  of  animal  herds;  detection  by 
several  sensors  simultaneously  is  necessary  for  reliable  tracking,  buc  individual 
sensors  have  a  limited  radius  of  sensitivity  (or  effective  horizon  distance),  and  may 
interfere  with  one  another  if  sited  too  close  together.  Since  trade-offs  may  be  possible 
becween  sensitivity,  interference  and  number  of  simultaneous  detections  necessary,  it 
would  be  desirable  to  examine  the  relationships  between  these  quantities  for  different 
types  of  cesselacion;  there  are  thus  direct  parallels  with  che  earth-station  and 
satellite-constellation  problems. 

We  first  consider  the  three  regular  tesselations,  as  shown  in  Fig  2.  The  basic 
dimension  for  each  is  che  edge-length  (or  inter-sensor  distance)  L  .  Each  tesselation 
may  be  divided  into  equal-area  cells  each  containing  a  single  sensor,  as  shown  on  che 
right-hand  part  of  each  diagram;  their  areas  correspond  to  rectangles  of  sides 
L  *  */3/2  L  for  the  triangular  tesselation,  L  *  L  for  the  square  tesselation,  and 

3/2  L  x  1/3/2  L  for  the  hexagonal  tesselation.  Hence  the  average  number  of  sensors  per 

2 

unit  area,  which  we  shall  denote  by  S  ,  is  1.15470/L  for  the  triangular  tesselation, 

2  2 

I .00000/L  for  the  square  tesselation,  and  0.76980/L  for  che  hexagonal  tesselation;  or, 

for  a  fixed  density  of  sensors  per  unit  area,  che  inter-sensor  distances  L  are  in  the 

ratio  of  1.07457  for  che  triangular  tesselation  to  1.00000  for  the  square  tesselation 

and  0.87738  for  the  hexagonal  tesselation. 

When  evaluating  satellite  patterns  over  a  spherical  Earth,  we  are  usually  comparing 
patterns  containing  che  same  total  number  of  satellites  (T) ,  distributed  over  a  fixed 
area  (the  spherical  surface  of  che  Earth),  on  che  basis  of  the  effective  horizon  dis¬ 
tance  (R^ax  n>  at  which  a  given  level  of  coverage  (n)  can  be  achieved,  and  possibly  also 

on  the  basis  of  the  associated  minimum  inter-satellite  distance  D  .  ;  R.  should 

Min  >lax,n 

preferably  be  small,  and  large.  We  may  apply  similar  concepts  to  the  tesselations, 

substituting  'sensors'  for  'satellites';  however,  since  there  are  no  potential  variations 
in  the  tesselations  corresponding  to  those  due  to  changing  satellite  phase  and  orbital 
inclination,  it  i3  sufficient  to  refer  only  to  values  of  Ry^  n  and  Dmj>;  • 

For  any  tesselation  T  and  n  will,  if  large,  be  closely  proportional  to  the 
associated  surface  area.  If  the  values  of  T,  n  and  the  associated  surface  area  are  all 
large  and  fixed  then  the  value  of  is  also  fixed,  independent  of  which  of  the 

tesselations  is  involved,  and  is  proportional  to  /n  .  On  che  other  hand  the  value  of  L 
(ta  of  Dydjj)  will  depend  on  which  tesselation  is  involved;  the  value  of  Dy^  will  be 
largest  (ie  most  favourable)  for  the  triangular  tesselation  and  smallest  for  the 
hexagonal . 


Tor  large  values  of  T  and  n  ,  circur.circ  les  defining  the  effective  radius  of 
coverage  will  enclose  large  numbers  of  the  unit  equal-area  cells,  so  that  there  will  be 
close  correspondence  between  the  measured  value  and  the  average  value  for  number  of 
sensors  enclosed  within  a  circle  of  given  radius,  and  hence  area.  However,  as  the  values 
of  I  and  n  become  smaller,  these  values  will  increasingly  tend  to  vary  from  the 
average  values  common  to  the  three  tesse lations ,  and  to  vary  in  different  ways  for  the 
three  tesse lations ,  dependent  on  the  detailed  geometry  of  each.  The  examination  of  this 
effect  is  discussed  in  the  following  section. 

3.2  Determination  of  R., _ „  for  tesselations 

_ >IAX ,  n _ _ 

We  consider  first  the  triangular  tesselation,  as  shown  in  Fig  3a.  For  this  case 
^MAX  1  corresponds  to  the  radius  of  the  circumcircle  of  three  adjacent  sensors  forming 
one  of  the  basic  triangles,  i*  the  circle  A;  its  centre,  also  marked  A,  lies  on  the 
central  vertical  line,  and  its  radius  may  be  calculated  to  be  0.577dL.  ^  corres- 

ponds  Co  che  radius  or  the  circumcircle  enclosing  one  other  sensor,  :'e  the  circle  B, 
with  radius  L;  this  actually  passes  through  a  total  of  six  sensors. 

To  develop  a  complete  list  of  values  of  ^  as  n  is  increased,  we  need  (i)  to 

identify  all  the  different  circles  which  pass  through  three  or  more  sensors;  (ii)  to  find 
for  each  circle  the  number  (x)  of  sensors  through  which  it  passes  and  che  number  (y)  it 
encloses,  thereby  obtaining  (as  in  section  2.2)  the  corresponding  value  or  values  of  n  ; 
(iii)  to  calculate  the  radius  of  each  circle;  and  then  (iv),  for  each  value  of  n  ,  to 
find  which  circle  has  the  largest  radius.  To  ensure  that  all  circles  have  been 
identified,  it  is  advisable  first  to  compile  a  list  of  all  which  pass  through  a  pair  of 
adjacent  sensors,  as  in  Fig  3a,  and  then  to  repeat  this  process  for  other  pairs  of 
sensors  situated  at  increasing  distances  apart;  one  example  for  another  pa'-  of  sensors 
is  shown  in  Fig  3b.  Each  list  may  be  checked  for  internal  consistency  in  terms  of  che 
increased  numbers  of  sensors  lying  on  or  within  successive  circles  of  increasing  radius, 
and  the  lists  may  then  be  combined  into  a  single  master  list  arranged  in  increasing  order 
of  radius,  the  duplications  which  inevitably  occur  between  different  lists  providing  a 
second  check  on  consistency.  (It  may  be  noted  that  circles  E,  G  and  -  are  duplicated 
as  between  Figs  3a  and  3b;  the  circle  identification  letters  were  allocated  after  compila¬ 
tion  of  the  master  list.)  Table  I  contains  as  examples  the  individual  lists  corresponding 
to  Figs  3a  and  3b,  and  Table  2  concains  the  master  list  for  the  triangular  tesselation;  in 
this  a  stepped  line  has  been  drawn  between  the  values  of  n  to  separate  the  critical  case 
'  that  corresponding  to  the  largest  circumcircle  radius)  for  each  value  of  n  ,  lying 
below  and  to  the  left  of  the  line,  from  non-critical  cases  lying  above  and  to  its  right. 

It  can  be  seen  from  Fig  2  that  the  hexagonal  tesselation  is  identical  to  the 
triangular  tesselation  with  every  third  sensor  in  each  row  removed.  Each  individual 
list  land  the  master  list)  therefore  contains  the  same  circle  radii  as  for  the  triangular 
tesselation,  but  with  different  values  of  x  and  v  ;  since  the  'holes'  will  appear 
in  different  places  for  different  alignments  of  the  tesselation,  some  circles  may  appear 
ir.  the  lists  with  more  than  one  pair  of  values  of  x  and  y  ,  while  others  may  not 
appear  at  all.  The  square  tesselation,  on  the  other  hand,  involves  a  different  set  of 
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circle  radii.  The  resulcing  master  lists  appear  in  Table  3  for  the  square  tesselation 
and  Table  4  for  the  hexagonal  tesselation. 

In  Tables  1  to  4  the  circumcircle  radii  have  been  expressed  as  a  function  of  L  , 
whereas  to  compare  them  in  the  manner  normally  used  for  satellite  constellations  they 
should  be  expressed  as  a  function  of  S  ,  che  average  number  of  sensors  per  unit  area, 

2  .  2 

which  we  found  to  equal  1.15470/L  for  the  triangular  tesselation,  1.00000/L  for  the 

2 

square  tesselation,  and  0.76980/L  for  the  hexagonal  tesselation.  For  a  circle  of 
radius  R  enclosing  a  large  number  of  unit  cells  of  a  tesselation  we  woul^  expect  the 


to  be  equal  to  iR  S  ,  ie  we  would  expect  R^  to  be 


corresponding  degree  of  coverage  n 

equal  Co  >'n/ itS  .  It  is  of  interest  to  compare  the  values  of  n  found  for  the 

three  regular  tesselations  for  relatively  small  values  of  n  with  these  expected  values. 

A  normalised  value  of  n  •  which  we  shall  denote  by  on  ,  can  be  obcained  by 

multiplying  the  values  of  radius/L  in  Tables  2  to  4  by  IVfS/n  ,  ie  by  1.07457»T7n 

for  the  triangular  tesselation,  by  /*/ n  for  the  square  tesselation,  and  by  0.87738/r/n 

for  the  hexagonal  tesselation:  values  of  o  would  be  expected  to  approach  1.0  for  large 

n 

,  and  hence  of  p  ,  should  preferably  be  small,  the 
JIAa,  n  n 

value  of  l/o  may  be  regarded  as  a  figure  of  merit  for  che  tesselation  in  respect  of 
n 

the  actual  value  of  R^^  n  achieved,  a  value  less  than  1.0  being  inferior  to  the 
'expected'  value. 


values  of  n  .  As  values  of 


Table  5  lists,  for  each  value  of 


the  critical  values  of  radius/L  for  the  three 


regular  cesselacions  as  found  in  Tables  2  Co  4,  with  the  corresponding  values  of  p  and 


I  /a 


In  Fig  4  values  of  1 


l/on  for  these  tesselations  have  been  plotted  against  n 


on  a  semi-logarithmic  scale;  this  suggests  that,  to  a  first  approximation,  values  o  f  I  /o^ 

for  these  tesselations  approach  the  value  1.0  exponentially  as  n  increases,  according 

to  che  relationship  1/p  •  I  •  pe  ,  where  p  -  0,16  and  q  »  0.08.  The  scatter  about 

the  straight  line  corresponding  to  this  relationship  appears  random  (though  in  reality 

all  individual  values  are  strictly  defined)  and  is  within  limits  of  factors  between  3.2 

and  1/3.2  on  the  value  of  p  ,  ie  between  p  »  0.05  and  0.512.  It  may  be  noted  that  if 

we  choose,  for  each  value  of  n  from  1  to  16,  whichever  of  che  three  tesselations 

provides  the  largest  value  of  1/p^  ,  the  resulting  values  (seven  from  the  triangular, 

seven  from  the  square  and  two  from  the  hexagonal  tesselations)  lie  becween  narrower 

limits  corresponding  to  p  »  0.05  and  p  »  0.17,  with  a  likely  value  of  p  of  about  0.092. 

(Mote  that  p  is  equal  to  the  value  of  1  -  l/o  when  n  ■  0.) 

n 

We  have  so  far  considered  only  the  regular  tesselations;  however,  their  very 
regularity  makes  them  appear  remote  from  the  r mditions  encountered  with  practical 
satellite  constellations,  which  have  numerous  different  values  of  inter-satellite  dis- 


D„.  must  be  found,  and  numerous  different  values  of 
Mm 


from  which 


must  emerge  for  any  given  value  of 


To  see  how  a  less  regular 


tance  from  which 

the  value  of  R,. 

>!ax,n 

pattern  affects  the  results  obtained,  three  other  tesselations  have  been  examined;  they 
were  chosen  to  show  the  effects  of  including  differing  face  shapes,  differing  vertex 
configurations  and  differing  edge-lengths.  Like  che  hexagonal  tesselation,  they  are  all 
obcained  from  che  triangular  tesselation  by  omitting  certain  sensors,  so  chat  the  same 
list  of  calculated  values  of  radius/L  mav  be  used. 


'll  I 


i  3 

The  first  of  these,  shown  in  Fig  5a,  we  shal  .  afer  to  as  'tri-hex  A’;  it  was  listed 
a 

by  Critchlow  as  seni-regular  tesselation  No.  I,  having  a  single  edge-length,  a  single 
vertex  configuration,  and  two  different  (triangular  and  hexagonal)  face  shapes.  It  is 
obtained  from  the  triangular  tesselation  by  omitting  no  sensors  from  one  row,  alternate 
sensors  from  the  next,  none  from  the  next,  and  alternate  sensors  from  the  following  row 
staggered  relative  to  those  in  the  second  row,  ie  omicting  one  sensor  in  four  overall. 

The  sensor  density  S  is  therefore  three-quarters  of  that  for  the  triangular  tesselation, 
is  0.86603/L",  and  the  relative  inter-sensor  distance  is  0.93060. 

The  second  of  these  other  tesselations,  shown  in  Fig  5b,  we  shall  refer  to  as 

9 

'tri-hex  B ’ ;  it  was  listed  by  Critchlow  as  demi-regular  tesselation  No. 6,  having  a 
single  edge-length,  two  different  vertex  configurations,  examples  of  which  are  ringed  on 
Fig  5b  (one  has  similar  faces  opposite,  the  other  similar  faces  adjacent),  and  two 
different  face  shapes,  again  triangular  and  hexagonal.  Like  tri-hex  A  it  is  obtained 
from  the  triangular  tesselation  by  omicting  alternate  sensors  from  even-numbered  rows, 
but  with  no  staggering  between  those  rows;  the  sensor  density  and  relative  inter-sensor 
distance  are  the  same  as  for  tri-hex  A. 

The  third  such  tesselation,  which  we  shall  refer  to  as  "the  rectangular  tesselation", 
is  shown  in  Fig  5c;  it  is  obtained  from  the  triangular  tesselation  by  omitting  alternate 
rows  entirely.  It  thus  has  two  different  edge  lengths,  L  and  L’  ,  with  L’  -  /3  L; 

L  ,  as  the  smaller  value,  corresponds  to  •  The  sensor  density  is  half  that  for  the 

triangular  tesselation,  i e  0.57735/L2,  and  the  relative  value  of  L  is  0.75984. 

As  with  the  hexagonal  tesselation,  the  'holes'  in  these  tesselations  appear  in 
different  places  for  different  alignments,  so  multiple  circumcircle  lists  are  required; 
Table  6,  corresponding  to  part  of  Table  1,  gives  as  an  example  the  lists  for  the  circum- 
circles  shown  in  Fig  3b.  Such  lists  are  then  combined  into  master  lists,  as  was  done  in 
Tables  2  to  4,  and  from  these  emerges  Table  7,  which  corresponds  to  Table  5  for  the 
regular  tesselations;  the  values  of  radius/L  in  Table  7  come  direct  from  the  master 
lists,  and  the  values  of  are  obtained  by  multiplying  the  values  of  radius/L  by 

0.93060*^/n  for  the  tri-hex  tesselations  and  by  0.75984»'-/n  for  the  rectangular 
tesselation. 

The  resulting  values  of  l/c  are  plotted  in  Fig  6;  for  comparison,  the  lines 

n 

corresponding  to  the  likely  and  limiting  values  for  the  regular  tesselations,  as  found 
from  Fig  4,  are  superimposed  on  this  figure.  Though  the  differences  are  not  marked,  it 
appears  that  the  likely  values  of  l/cn  for  the  tri-hex  tesselations  may  be  slightly  lower 
chan  chose  for  the  regular  tesselations  for  small  values  of  n  (though  they  would  still 
approach  1.0  for  large  values  of  n  );  if  so,  this  may  be  associated  with  the  larger  size 
of  the  unit  equal-area  cells  of  the  tri-hex  constellations,  each  enclosing  three  sensors, 
as  shown  on  the  right-hand  side  of  Fig  5a&b.  More  nociceablv,  it  appears  that  the  likely 
values  of  I/c^  for  the  rectangular  tesselation  may  also  be  lower;  this  may  be  associated 
with  the  fact  that  the  value  of  L  is  now  less  than  the  mean  inter-sensor  distance.  In 
the  relationship  l/c  »  I  -  pe  qn,  the  value  q  -  0.08  still  appears  appropriate,  but  the 
likely  value  of  p  for  these  tesselations  appears  to  be  about  0.21,  with  limiting  values 
of  p  ■  0.21  x  2.4  and  0.21/2.4,  0.0875  and  0.504. 


A  point  of  interest  to  note,  for  the  three  regular  tesselations  in  Table  5  and  the 
tri-hex  tesselations  in  Table  7,  is  that  when  n  is  2  the  value  of  /L  is  *-0 

in  each  case,  is  R^^  ,  -  L;  this  is  because,  with  a  single  value  of  edge-length,  the 
(n  »  2)  circumcircle  must  have  a  vertex  at  its  centre,  and  its  radius  is  therefore  equal 
to  the  edge-length.  Hence  DMIN  «  For  the  rectangular  tesselation,  however, 

with  its  differing  edge-lengths,  is  less  chan  2  ■ 

3.3  Determination  of  R^^,  for  regular  polyhedra 

Having  established  some  of  the  trends  associated  with  patterns  based  on  regular 
tesselations  of  a  plane  surface,  it  is  worth  considering  the  regular  polyhedra  as 
representing  a  basis  for  patterns  covering  a  spherical  surface,  and  hence  as  an  inter¬ 
mediate  step  between  the  tesselations  and  practical  constellations  of  orbiting 
sacellices. 

The  vertices  of  the  regular  polyhedra  may  be  considered  as  points  on  a  spherical 
surface  at  latitudes  and  longitudes  (for  one  possible  orientation)  as  follows: 

Tetrahedron  :  latitude  35.2644°  ac  longitudes  0°  and  180°; 

latitude  -35.2644°  at  longitudes  90°  and  270°. 

Occahedron  :  latitude  35.2644°  at  longitudes  0°,  120°  and  240  ; 

latitude  -35.2644°  at  longitudes  60°,  180°  and  300°. 

Cube  :  latitudes  35.2644°  and  -35.2644° 

at  longitudes  0°,  90°,  180  and  270  . 

Icosahedron  :  latitudes  52.6226°  and  -10.8123°  at  longitudes  0°,  120°  and  240°; 

latitudes  10.8123°  and  -52.6226°  at  longitudes  60°,  180°  and  300°. 

Dodecahedron  :  latitudes  52.6226°  and  10.8123°  at  longitudes  0°,  72°,  144°,  216° 

and  288°;  latitudes  -10.8123°  and  -52.6226°  at  longitudes  36°,  108°, 
180°,  252°  and  324°. 

From  these  coordinates  one  may  (using  standard  formulae  of  spherical  trigonometry) 
calculate  circumcircle  radii,  in  terms  of  the  angle  subtended  at  the  centre  of  the 
Earth,  and  for  each  circle  find  the  corresponding  values  of  x  and  y  ,  and  hence  the 
relevant  values  of  n  .  These  are  listed  in  Table  8;  brackets  have  been  placed  round 
non-critical  values  of  n  .  As  for  the  regular  tesselations,  L  (or  D  jjj)  is  equal  to 
the  second  radius  listed  in  each  case. 

Values  of  Rj^  n  extracted  from  Table  8  are  listed  in  Table  9  for  all  values  of 
n  relevant  to  each  of  che  regular  polyhedra.  The  procedure  for  deriving  values  of  on 

is  a  little  different  from  the  case  of  the  tesselations.  Here  we  have  a  fixed  area, 

2 

that  of  che  spherical  surface  (4tr  ),  associated  with  the  total  number  (T)  of  vertices  of 

1  _ _ 

the  polyhedron,  so  that  S  *  T/drr".  Since  0  *  k-rS/n  <  R,,.„  ,  values  of  c  for  the 

n  MAa  ,  n  _  n 

different  polyhedra  are  obtained  by  multiplying  R^^  n  ''/T/2r*n  ,  where  2r 

corresponds  to  360°/f,  is  by 
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tor  the  tetrahedron 


:  5 


180 
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for  the  octahedron, 
for  the  cube, 

for  the  icosahedron  and 

for  the  dodecahedron. 


The  resulting  values  of  c  and  I /■  are  listed  in  Table  9,  and  values  of  \/z  are 
s  n  n  n 

plotted  against  n  in  Fig  7a,  with  the  straight  lines  from  Fig  4  corresponding  to  likely 
and  limiting  values  for  regular  tesselations  superimposed  for  comparison.  While  all  the 
calculated  values  of  1 /on  are  plotted  on  Fig  7a,  we  are  really  only  interested  in  those 
derived  from  values  of  less  than  90°,  since  larger  values  give  no  visibility 

round  the  spherical  Earth;  these  lower  values  are  therefore  distinguished  by  solid 
symbols.  It  appears  from  Fig  7a  that  the  likely  values  of  l/o^  £°r  regular  polyhedra 
may  well  be  following  a  similar  trend  to  that  for  regular  tesselations  (chough  at 
slightly  lower  values)  while  [kIAX  n  as  less  than  90°,  but  that  as  n  exceeds  90° 

and  approaches  180°  the  values  of  1/c^  rend  to  fall.  While  the  usual  scatter  of  values 
makes  it  difficult  to  identify  such  trends  with  confidence,  it  appears  that  the  relation¬ 
ship  l/on  «  1  -  pe  qn  may  again  be  appropriate  when  n  is  less  than  90°,  with  q 

again  equal  to  0.08  and  the  likely  value  of  p  equal  to  0.21,  with  limiting  values  of 
p  »  0 . 2 1  *  1.8  and  0.21/1.8. 

3.4  Determination  of  for  tesselations  and  polyhedra 

We  have  already  noted  values  of  L  ,  which  also  represent  the  values  of  > 

for  the  various  tesselations  and  the  regular  polyhedra;  L  is  proportional  to  1 / /s  , 
and  relative  values  are,  for  the  tesselations,  1.07457  for  the  triangular,  1.00000  for 
Che  square,  0.87738  for  Che  hexagonal,  C. 93060  for  tri-hex  A  and  B  and  0.75984  for  the 
reccangular;  and  for  the  regular  polyhedra,  1.07796  for  the  tetrahedron,  1.08540  for  the 
octahedron,  0.98216  for  the  cube,  1.08191  for  the  icosahedron  and  0.92060  for  the 
dodecahedron.  Thus  the  four  triangular-faced  figures  all  provide  values  approximating 
to  1.1,  compared  with  1.0  for  the  two  square-faced  figures,  0.92  for  the  pentagonal-faced 
dodecahedron  and  0.88  for  che  hexagonal  tesselation.  The  rectangular  cesselation,  with 
two  different  edge- lengchs ,  gives  a  considerably  lower  value  at  0.76,  while  tri-hex  A 
and  B  fall  between  the  levels  for  triangles  and  hexagons. 


It  appears  appropriate  to  normalise  che  values  of  D,^,,  >  as  was  done  for  R>^y  n 
so  as  to  give  an  expected  value  of  1.0  for  a  large  area  of  an  optimum  tesselation,  and 


smaller  values  for  less  satisfactory  values  of  D, 


MIN 


Since  should  preferably  be 


large,  and  the  triangular  tesselation  is  the  one  giving  the  largest  value,  normalised 


16 


values  of  D, T„  (which  we  shall  call  A  )  raav  be  obtained  simplv  by  dividing  the 
MIN 

relative  values  by  1.07457.  Resulting  values  of  1  are: 

triangular  tesselation  1.0000,  tetrahedron  1.0032,  octahedron  I .0101 , 
icosahedron  1.0068; 

square  tesselation  0.9306,  cube  0.9140; 
dodecahedron  0.8567; 
hexagonal  tesselation  0.8165; 
tri-hex  A  and  B  0.8660;  and 
rectangular  tesselation  0.7071. 

The  fact  thac  the  triangular-faced  polyhedra  have  values  of  A  slightly  exceeding  1.0 
may  detract  from  the  neatness  of  this  approach,  but  is  irrelevant  to  our  ultimate 
objective  of  deriving  values  of  ant*  *  for  practical  satellite  constellations. 

Values  of  X  may  be  obtained  directly  from  values  of  which  have  been 

measured  on  a  spherical  surface  in  terms  of  degrees  subtended  at  the  centre  of  the  sphere 
by  multiplying  by  At/37  ij  360°  »  /t/218.254. 

3.5  A  combined  criterion 

In  the  light  of  our  consideration  so  far  of  the  application  of  the  separate 
criteria  n  and  ,  and  the  corresponding  normalised  criteria  1/;^  and  '•  , 

it  is  worth  examining  the  possibility  of  a  combined  criterion  to  meet  the  objective 
suggested  in  section  2.4. 

An  obvious  candidate  for  consideration  is  the  product  of  the  two  normalised 

criteria  '/°n  and  X  ;  since  each  approaches  the  value  1.0  under  ideal  conditions, 

while  otherwise  giving  smaller  values,  their  product  will  do  the  same.  '■/;  is 

n 

proportional  to  D  ../ft,  ;  this  is  non-dimensional,  since  both  DWT„  and  R,,,„ 

M  l N  *tAA  ,n  MIN  MAX  f  n 

are  distances  measured  on  the  surface  of  a  sphere.  We  shall  denote  the  product  « / o  by 
en  ;  since  ' /°n  "  360°  x  Zn/T/eR ,  and  \  «  *  At.,3/2/360o,  we  have 


„  —  M  [  \ 

-  0 • 52504f n  j  —  , 

MAX,n 

and  we  shall  test  the  use  of  this  prospective  criterion  when  we  come  co  examine  practical 
satellite  constellations.  For  the  tesselations  and  regular  polyhedra,  we  need  merely 
note  that  combining  1/c^  and  ■  in  this  manner  serves  to  distinguish  between  the  three 
regular  tesselations,  giving  the  advantage  to  che  triangular  tesselation,  and  that  it 
likewise  gives  the  advantage  to  the  three  triangular-faced  regular  polyhedra  over  che 
cube  and  the  dodecahedron. 


APPLICATION'  OF  ANALYSIS  TO  DELTA  PATTERNS 


4.  1  Derivation  of  K,.v  and  D. , . 

_ MAX ,  n _ Min 

In  Ref  3  the  author  listed  satellite  patterns,  containing  a  total  number  of 
satellites  (T)  between  5  and  15,  which  for  each  value  of  T  appeared  most  economically 
to  provide  single,  double,  triple  or  quadruple  coverage  of  the  whole  surface  of  the 
Earth.  These  patterns  all  belonged  to  a  class  of  pattern  which  the  author  had  pre¬ 
viously'  described  as  'delta  patterns'  (and  which  have  since  been  referred  to  elsewhere'^ 
as  ’Walker  constellations').  As  noted  in  section  I,  a  delta  pattern  consists  of  equal 
numbers  of  evenly-spaced  satellites  in  each  of  P  evenly-spaced  and  equally-staggered 
planes,  all  having  the  same  inclination  5  to  a  reference  plane  (usually,  but  not 
necessarily,  the  equator).  It  may  be  identified  by  the  three-integer  code  reference 
T/P/F,  where  F  (measured  in  units  of  360°/T)  is  a  measure  of  the  relative  phasing 
between  satellites  in  adjacent  planes;  specifically,  when  a  satellite  in  one  plane  is  at 
its  ascending  node,  a  satellite  in  the  adjacent  plane  having  a  more  easterly  ascending 
node  is  (360  F/T)°  past  its  ascending  node. 


Table  10  lists,  for  each  value  of  T  ,  the  delta  pattern  giving  the  smallest  value 

of  n  (for  values  of  n  from  I  to  A),  as  quoted  in  Ref  3,  together  with  this  value 

of  RMAX  n  ,  the  inclination  J  at  which  it  occurred,  and  the  corresponding  value  of 

D.,.  From  these  values  the  normalised  criteria  c  ,  l/o  ,  '  and  e  (defined  as  for 
Mn  n  n’  n 

the  regular  polyhedral  have  been  calculated,  and  are  also  listed. 


It  should  be  noted  that,  for  T  «  10,  the  pattern  10/10/7  actually  gives  a  slightly 
smaller  value  of  RJ1AX  |  (51.5°  at 


5  of  47.9  ,  with  D  ■  0  )  than  the  pattern 
opt  min. 


10/5/2  listed  here.  Pattern  10/10/7  has  been  listed  by  Mozhaev1'  and  Ballard'"",  but  the 
present  writer  chose  to  exclude  from  his  listings  any  pattern  giving  a  value  of  D^n 
less  than  3°,  as  unlikely  to  be  acceptable  for  a  practical  satellite  system;  10/10/7  is 
a  patcern  for  which  P  and  F  -  T/P  are  both  even  numbers,  which  identifies  it  (see 


section  3.5  of  Ref  4)  as  one  for  which  D, 


Min 


0  at  all  inclinations,  and  so  pattern 


10/5/2  appears  preferable.  Pattern  1 0/ 10/7  would  have  zero  values  of  1  and  e 


Table  10  also  includes  five  further  patterns,  having  values  of  T  between  16  and 
24,  which  were  listed  in  Ref  4  as  giving  optimum  values  of  R^^  n  for  values  of  n 
between  4  and  7. 


Fig  7b  is  a  semi-log  plot  against  n  of  the  values  of  l /c ^  corresponding  to 
these  smallest  values  of  n  ,  again  with  the  lines  from  Fig  4  (corresponding  to 

likely  and  limiting  values  for  regular  tesselations)  superimposed  on  it.  It  appears  that 
the  value  of  q  -  0.08  is  again  appropriate,  with  limiting  values  of  p  of  about  0.20 
and  0.34,  and  a  likely  value  of  p  of  about  0.26.  Figs  7c  and  "d  will  be  discussed  in 
the  following  section. 

-.2  Derivation  of  D,,.,.  and  R 

_ MIN _ Max.n 

Prior  to  197'  the  author  did  not  attempt  to  determine  values  of  D  for  delta 
patterns,  concentrating  instead  on  determination  of  values  of  R^  ..  and  c^e  associated 


values  of  D  .  .  In  1977,  however,  following  up  an  observation  noted  in  Ref  4  that  two 

Mm 

series  of  patterns  listed  in  its  Table  1  appeared  to  give  the  largest  values  of 
for  values  of  T  above  and  below  10  respectively,  values  of  and  the  inclinations 

(iQpt  q)  at  which  they  occurred  were  determined  (using  an  abbreviated  version  of  the 
program  COCO)  for  four  such  series  of  patterns,  for  values  of  T  up  to  25,  and  a  plot  of 
^MIN  a8aanst;  T  was  provided  in  Fig  3  of  a  brief  paper^  presented  in  1978.  This  i's 
reproduced  as  Fig  8  of  the  present  Report,  while  in  Fig  9  values  of  Q  are  plotted 

against  T  for  the  same  four  series  of  patterns.  These  patterns,  identified  by  their 
code  reference  T/P/F,  are  listed  in  Table  11  together  with  the  corresponding  values  of 

°MIN*  5opt,D  and  X  • 

The  four  series  of  patterns  are  those  which  would  follow  single  non-self-crossing 
Earch-tracks  if  used  at  periods  of  12  h,  16  h,  18  h  or  19.2  h,  ie  they  would  complete  L 
orbits  in  M  days,  where  L:M  equals  2:1,  3:2,  4:3  or  5:4.  It  was  shown  in  Ref  4  that 
non-self-crossing  Earth-tracks,  corresponding  to  a  well-spaced  distribution  of  satellites, 
are  obtained  when  L  -  M  »  1 ,  and  that  this  favourable  distribution  is  retained  independent 
of  the  period  at  which  the  pattern  is  actually  used.  The  delta  patterns  associated  with 
these  values  of  t:M  are  identified  by  having  P  equal  to  the  value  of  T  divided  by 
the  highest  common  factor  of  T  and  M  ,  and  F  equal  to  T(kP  -  L)/PM  ,  where  k  takes 
whatever  integer  value  is  necessary  to  make  F  an  integer  in  the  range  0  to  P  -  I . 


The  largest  value  of  is  given  by  the  2:  1  series  for  values  of  T  from  5  to 

9,  by  che  3:2  series  from  10  to  24,  and  by  Che  4:3  series  from  25  upwards;  the  second- 
largest  value  is  given  by  che  3:2  series  for  values  of  T  below  10  and  above  24,  by  the 
2:1  series  from  10  to  15  and  by  the  4:3  series  from  16  to  24.  Fig  8  suggests  that  the 
5:4  series  would  become  relevant  when  T  reached  values  in  the  upper  twenties  or  lower 
thirties.  These  four  series  of  patterns  generally  provide  the  largest  values  of  D 

min 

for  all  the  patterns  having  a  particular  value  of  T  ;  for  example,  a  check  of  all  delta 
patterns  having  T  »  13  shows  that,  while  these  four  have  values  of  of  44.1°, 

47.3°,  41.7°  and  35.6°,  che  pattern  in  the  8:5  series  has  a  value  of  31.7°  and  the 
remainder  all  have  values  below  30°.  For  small  values  of  T  ,  some  patterns  reach  their 
largest  value  of  at  zero  inclination  (at  which  R^^  n  *  90°);  these  values  all  lie 

on  the  broken  line  on  the  left  side  of  Fig  8.  Otherwise,  cross-reference  to  Fig  9  shows 
that  the  largest  values  of  occur  at  inclinations  between  51°  and  66°.  Reference 

to  Table  10  shows  that  there  is  not  the  same  regular  pattern  to  che  occurrence  of  minimum 
values  of  R^^^  n  •  though  about  802  of  the  values  of  t  listed  fall  within  the  same 
range  of  inclinations. 


As  an  example  of  the  conditions  under  which  occurs  for  che  patterns  giving 

non-self-crossing  Earth-tracks,  Fig  10  shows  pattern  13/1 3/5  at  the  inclination  of  57.6 

corresponding  to  ^  and  at  phase  :  ■  1.0;  under  these  conditions  che  distances 

between  satellites  CH  and  LN  are  both  equal  to  DUT,,  .  If  5  were  increased  then 

MIN 

D„.  (associated  with  CH)  would  be  less,  while  if  5  were  reduced  then  D„. 

Min  Min 

(associated  with  LN)  would  be  less.  For  ocher  such  patterns  may  occur  for  both 

pairs  when  :  ■  0  or  for  one  pair  when  s  »  0  and  for  the  other  when  :  •  1.0;  buc  in  all 
such  cases,  one  pair  is  symmetrically  placed  near  the  closed  end  of  a  loop  and  the  other 
near  the  open  end  of  a  loop. 


and 


were  not  extended 


The  i 9 7 7  computer  runs  to  determine  values  of  D 


to  determine  corresponding  values  of  R,, 

>lax ,  n 


MIN  ' op  t ,  D 

However,  in  many  cases  it  has  been 
possible  to  refer  back  to  the  printouts  from  the  runs  of  the  program  COCO  which  deter¬ 
mined  values  of  R.,,..  for  pacterns  with  T  between  5  and  15  in  1971-2,  and  chose  for 
MAX ,  n 

some  pacterns  with  T  between  16  and  25  in  1974-5,  and  by  interpolation  obtain  a  good 

approximation  to  the  values  of  R,,  corresponding  co  inclination  5  .  .  This  has 

Max,n  opc,D 

made  it  possible  to  estimate  values  of  c  ,  1/c  and  £  ,  as  well  as  a  ,  corresponding 

r  n  n  n  ’ 

to  that  inclination.  Values  are  listed  in  Table  12  for  the  patterns  providing  the 
largest  values  of  D  ^  for  eac^  value  of  T  ,  and  in  Table  13  for  those  providing  the 
second-largest  values. 

While  the  printouts  were  being  examined  and  values  of  £  calculated  at  £  and 

n  opt 

£  _  ,  values  were  also  calculated  at  other  inclinations  in  order  to  find  the  maximum 

opt,D  ’ 

values  of  £  and  the  inclinations  at  which  they  occurred.  It  was  found  that  (with  one 


exception)  the  maximum  value  of  c  alwavs  occurred  at  £  _ 

e  '  n  opt,D 


This  is  because  D,.. 

Min 


always  falls  off  much  more  rapidly  as  £  is  either  increased  or  decreased  from  the  peak 

at  i  „  than  R^,  increases  as  £  is  either  increased  or  decreased  from  * 

opc,D  Max,n  opt 

(as  mav  be  seen,  for  example,  from  Figs  6  and  7  of  Ref  2);  the  value  of  e  is  thus 
'  n 

much  more  sensitive  co  0M.  chan  co  R^,  .  The  single  exception  was  patcern  5/5/1, 

Min  .  .sx ,  n 

which  lies  on  che  broken  line  in  Fig  8.  For  this  pattern  occurs  at  zero  inclina¬ 
tion  and  the  value  of  decreases  only  slowly  as  £  is  increased;  the  maximum  value 


therefore  occurs  at 


,  as  shown  in  Table  10. 


n  opt 

The  values  of  in  Fig  8,  excluding  those  lying  on  the  broken  line,  are 

replocted  in  Fig  1 1  in  the  form  of  a  log-log  plot  of  X  against  T  (with  expanded  X 
scale).  This  shows  that,  as  T  is  increased,  values  of  X  corresponding  to  each  of  che 
(L  -  M  «  I)  series  tend  to  rise  co  a  peak  and  then  fall  away  again,  crossing  as  they  do 
so  the  rising  values  for  each  of  che  subsequent  series  as  these  in  turn  approach  their 
own  peaks.  A  line  representing  the  envelope  of  the  peaks  and  another  representing  the 
lower  envelope  of  che  intersections  of  the  curves  representing  consecutive  series, 
corresponding  to  che  broken  lines  in  Fig  II,  indicate  the  upper  and  lower  limits  of  the 
largest  values  of  «  for  a  given  value  of  T  ,  and  over  the  range  of  values  of  T  for 
which  results  are  available  ic  appears  Chat  these  can  be  adequately  represented  on  this 
log-log  plot  by  che  straight  lines  corresponding  to  •  •  1.01  T  and  •  «  0.96  T 
with  a  likely  value  between  these  limits  represented  by  the  straight  line  X  ■  0.99  T 


For  che  delta  pacterns  from  che  (L  -  M  *  I)  series  giving  che  largest  values  of 

D,,tv  (and  hence  of  ■  )  for  each  value  of  T  ,  the  values  of  1 !:  at  £  _  are 

..IN  n  opt.D 

plotted  against  n  in  Fig  7c;  for  chose  patterns  giving  che  second-largest  values  of 
0,,^,.  ,  values  of  l/c  are  plotted  in  Fig  7d.  As  in  Fig  ~b  it  appears  that,  over  the 
limited  range  of  values  of  n  for  which  results  are  available,  the  value  of  q  •  0.08 
remains  acceptable;  for  Fig  7c  it  appears  that  the  likely  value  of  p  is  about  0.32,  with 
limits  of  0.32,1.33  -  0.24  and  0.32  *  1.33  ■  0.-3,  while  for  Fig  7 d  che  likely  value  is  a 
liccle  higher  and  che  limits  somewhat  wider.  Thus  optimising  for  rather  than  for 

R,,^.  _  appears  cc  increase  che  value  of  p  such  chat  the  likely  value  approaches  what 
was  trie  less  satisfactory  of  the  limiting  values  for  R,^  n  ,  though  the  corresponding 
values  of  are  naturally  much  improved. 


4.3  The  combined  cricerion 


Values  of  e  ac  che  inclination  5  have  been  listed  in  Table  10,  and  values 
n  opt 

ac  i  „  in  Tables  12  and  13  for  the  relevant  patterns.  The  overall  average  of  the 
opt  ,D 

values  in  Table  12  is  0.566,  of  those  in  Table  13  is  0.529,  and  of  those  in  Table  10  is 

0.353.  It  has  already  been  noted  that  it  was  found  that,  in  all  cases  where  5  _ 

o  °PC>D 

occurred  at  an  inclination  acceptable  for  whole-Earch  coverage  and  not  at  0°,  the 

maximum  value  of  £  occurred  at  a  _  . 

n  opt.D 

When  it  was  decided  to  test  £fl  as  a  combined  criterion  for  use  with  practical 

satellite  patterns,  it  was  thoughc  that  ic  might  point  to  some  compromise  choice  as 

becween  those  patterns  which  gave  the  best  values  of  n  an<*  Chose  which  gave  the 

best  values  of  D„_„  ,  or  between  the  inclinations  S  and  6  -  for  a  single 

MIN  ’  opt  opt ,D  ° 

paccern.  Instead  the  test  has  shown  that  there  is  no  such  compromise;  che  choice  of 

pattern,  and  of  inclination  for  that  pattern,  is  so  much  more  sensitive  to  D  than 

MIN 

co  R,,.„  that,  if  it  is  accepted  that  D,,.,,  has  any  relevance  to  che  choice  to  be 
MAA ,  n  MIN 

made,  it  almost  inevitably  becomes  che  dominant  factor  in  thac  choice.  In  that  respect, 
the  exercise  may  be  considered  co  have  failed;  however,  the  information  it  has  provided 
regarding  the  importance  of  DMIN  is  of  considerable  significance.  is  very 

much  simpler  to  calculate  than  n  ,  and  if  it  is  possible  to  use  ic  as  the  sole 

criterion  during  the  p'reliminary  screening  of  satellite  patterns  for  a  particular 
application  then  much  cime  and  effort  can  be  saved. 

While  these  results  have  suggested  that  D,^  may  prove  a  more  generally  useful 

criterion  chan  n  .  ic  was  felt  worthwhile  to  examine  separately  the  values  of 

»,  l/on  and  cn  for  those  delta  patterns  giving  both  the  largest  and  the  second- 

largest  values  of  DUT„  for  each  value  of  T  ,  as  has  been  done  in  Tables  12  and  13. 

One  reason  for  doine  so  is  that,  near  che  cross-over  between  successive  (L  -  M  »  1) 

series,  che  values  of  A  will  be  very  similar  and  so  the  values  of  !/on  ®ay  determine 

which  pattern  has  the  larger  value  of  £  .  Another  reason  is  chat  the  advantages  of  a 

n 

pa'ticular  value  of  6  or  P  may  become  a  factor  in  an  otherwise  closely-balanced 
choice;  values  of  6opt  Q  differ  significantly  between  the  two  cables  (as  shown  also 

by  Fig  9),  and  for  most  values  of  T  divisible  by  2  or  3  the  value  of  P  differs 

between  che  patterns  in  these  two  tables. 

Comparing  results  for  similar  values  of  T  and  n  ,  ic  is  noted  that  when  T  »  24 

the  values  of  and  A  for  the  patterns  in  Tables  12  and  13  are  virtually  identical, 

though  those  for  24/12/9  are  actually  marginally  greater  than  those  for  24/8/4.  Their 

overall  similarity  is  emphasised  by  the  fact  that  the  values  of  1 /:  ,  and  hence  of  c  , 

n  n 

are  greater  for  24/12/9  when  n  1  I,  2  or  5,  but  are  greater  for  24/8/4  when  n  ■  3,  4,  6 

or  7.  Elsewhere  in  the  tables  there  are  only  two  cases  (T  *  1 1 ,  n  *  2  and  T  *  12,  n  *  I) 

where  a  larger  value  of  1 / o  results  in  the  value  of  c  being  greater  in  Table  13,  and 

n  n 

six  cases  (T  *  7,  n  *  2;  T  *  9,  n  *  3:  T  ■  12,  n  *  2  and  3;  T  *  13,  n  »  2;  and  T  *  14, 
n  *  2)  where  the  value  of  l/p^  is  larger  in  Table  13,  but  not  sufficiently  so  to  outweigh 
the  larger  value  of  \  in  Table  12;  in  most  cases  Table  12  contains  the  larger  values  of 
-n  as  well  as  of  X  .  We  would  therefore  expect  that  the  pattern  in  the  series  giving 
the  largest  value  of  should  be  given  first  consideration,  though  the  pattern 


hi  I 
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jiving  che  second- large  sc  value  might  be  a  suitable  choice  if  chere  were  other  reasons 
for  preferring  it. 

- • -  Discussion 

Reviewing  che  macerial  presenced  in  che  preceding  sections,  in  conjunction  with 
that  previously  presented  in  Ref  A,  suggests  the  following  overall  picture  of  the 
provision  of  multiple  whole-Earch  coverage  by  means  of  constellations  of  nigh-altitude 
circular-orbit  satellites  forming  delta  patterns. 

For  any  given  value  of  T  ,  che  total  number  of  satellites  in  the  pattern,  there 
will  be  a  substantial  number  of  different  delta  patterns  available  for  consideration; 
this  number  is  equal  to  the  sum  of  all  the  factors  of  T  ,  including  1  and  T  .  A 
proportion  of  these  patterns  will  give  values  of  R^^  n  which  may  be  considered  'good', 
while  a  further  proportion  will  give  values  which  may  be  considered  'poor';  the  methods 
described  in  section  3  of  Ref  U  will  identify  some  of  those  patterns  giving  'poor'  values 
of  n  and/or  of  •  Among  those  giving  'good'  values  of  R^^  n  •  3  similar 

situation  is  likely  to  exist  to  chat  illustrated  in  Fig  U  for  the  regular  tesselations, 
in  Fig  6  for  the  other  cesselations  and  in  Fig  7a  for  the  regular  polvhedra;  while  all 
remain  within  a  band  of  values  which  may  be  considered  'good',  there  is  an  irregular 
variation  with  the  value  of  n  of  their  relacive  order  of  merit,  determining  which  may 
be  considered  'best'.  The  particular  pattern  giving  the  smallest  value  of  R^^  n  for 
any  particular  values  of  T  and  n  can  be  found  through  a  detailed  investigation,  but 
is  otherwise  unpredictable;  while  the  'best'  values  of  n  fall  within  a  more 

narrowly  defined  band  than  the  'good1  values. 

Turning  temporarily  from  R^^  n  to  >  it  was  shown  in  Ref  5  (and  in  Fig  8  of 

che  present  Report)  chat  che  patterns  providing  the  largest  values  of  were  chose 

which  could  produce  single  non-self-crossing  Earth-tracks,  occurring  in  series  for  which 
L:M  is  such  that  L  -  M  ■  I.  Patterns  with  these  charac teris tics  usually  avoid  the 
adverse  features  which  lead  to  a  placing  in  the  'poor'  category  for  R^^  n  .  and 
instead  are  to  be  found  in  the  'good'  category,  with  a  fair  chance  of  appearing  as  'best' 
at  any  particular  value  of  n  .  Since  these  patterns  are  readily  identifiable  by  means 
of  formulae  developed  in  Ref  U  (and  reproduced  in  che  following  section) ,  they  provide  a 
convenient  starting-point  for  any  coverage  study;  moreover,  we  have  already  shown  that 
thev  provide  the  largest  values  of  our  test  criterion  z 

n 

While  these  conclusions  regarding  the  desirability  of  identifying  those  natcems 
providing  the  largest  values  of  D  j,  (or  '  )  have  been  developed  specifically  through 
examination  of  the  characteristics  of  delta  patterns,  it  might  be  expecced  that  they 
would  have  some  relevance  also  to  ocher  types  of  pactern  which  do  not  share  all  the 
characceris tics  of  delta  patterns. 

- .  5  Predictions 

The  results  wnich  have  been  discussed  in  che  preceding  sections  provide  new  means 
for  predicting  the  values  of  R,  ,  or  of  D,  ...  and  the  associated  value  of  Rv  , 

which  might  be  available  for  values  of  T  and  n  greater  than  those  which  have  already 


been  explored.  For  3.  modest  extrapolation,  and  with  somewhat  less  confidence  for  a 
larger  extrapolation,  we  may  use  the  predictions  that: 


W  <at  soPt>  *  °n  de8r 


where  I  /o^  »  I  -  pe  and  p  lies  between  limiting  values  of  0.20  and  0.34,  with 

an  expected  value  of  0.26;  hence 


n(i^) 


2  o' 


MAX,n 


Also 


D  (at  6  )  -  2-8‘i— -  degrees, 
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where  X  lies  between  limits  of  1.01  T  and  0.96  T  with  an  expected  value  of 

0.99  T0'1.  R,,  (at  4  _)  is  given  by  the  same  formulae  as  ,  but  with  p 

Max,n  op  c ,  JJ  7u\A,n 

lying  between  limiting  values  of  0.24  and  0.43,  with  an  expected  value  of  0.32. 

For  convenience,  values  of  l/o  for  the  six  relevant  values  of  p  are  presented 

2  n 

in  Table  14,  and  values  of  o  in  Table  15,  for  values  of  n  from  1  to  10;  also  the 
’  n 

three  values  of  D^N  are  presented  in  Table  16  for  values  of  T  from  26  to  40. 

No  means  are  available  of  predicting  either  the  value  of  associated  with  the 

value  of  R^^  n  at  ®opt  ’  °r  the  Parti-cuI-ar  pattern  which  would  provide  that  smallest 
value  of  R^^  n  •  It  can  be  predicted  that  che  largest  value  of  D».fXN  be  provided 

by  a  pactem  in  one  of  che  (L  -  M  •  1)  series  (identified  by  L:M)  ,  and  a  very  rough 

extrapolation  of  Fig  11  suggests  that  from  T  ■  25  to  about  30  it  would  be  from  the  4:3 
series  (or  perhaps  the  3:2  series);  from  T  ■  about  30  to  about  36,  from  the  4:3  series 
(or  perhaps  the  5:4  series);  from  T  -  about  36  to  about  42,  from  the  5:4  series  (or 

perhaps  che  4:3  series);  and  from  T  »  about  42  to  about  48,  the  5:4  series  (or  perhaps 

che  6:5  series).  The  specific  pattern  T/P/F  which  is  the  member  of  each  series  for  a 
particular  value  of  T  may  be  determined  from  che  following  formulae,  developed  in  Ref  4: 
P  is  equal  Co  T  divided  by  Che  highest  common  factor  of  T  and  M  ,  and  F  is  equal 

to  T(kP  -  L)/PM  ,  where  k  takes  whatever  integer  value  is  necessary  to  make  F  an 
integer  in  the  range  from  0  to  P  -  1 .  The  patterns  forming  the  3:2  series,  the  4:3 
series,  the  5:4  series  and  the  6:5  series  for  values  of  T  from  26  to  40  are  listed  in 
Table  17;  for  values  of  T  below  26  they  were  listed  in  Table  1  of  Ref  4  (and,  for  the 
first  three  of  these  series,  in  Table  II  of  che  present  Report). 

It  should  be  particularly  noted  that  4  -  is  the  inclination  ac  which  the  value 

opt.D 

of  D,^,.  largest,  but  that  varying  the  inclination  is  likely  to  reveal  a  smaller  value 
of  R\,ax  n  ■  Thus  if  the  pactern  which  provides  the  largest  value  of  at  a  parti¬ 
cular  value  of  T  does  not  meet  the  coverage  requirements  at  :  ^  g  because  it  has  coo 


large  a  value  of  n  at  chat  inclination,  it  is  still  worth  finding  the  value  of 

5,  for  that  pattern  at  :  to  see  if  the  requirements  can  then  be  met,  even  though 

the  value  of  is  reduced.  For  some  values  of  X  and  n  the  same  pattern  provides 

botn  Che  largesc  value  of  D,,.,.  at  i  _  and  the  smallest  value  of  R.,.,,  at  c  , 
s  MIN  opt,D  MAX.n  opt 

as  shown  by  comparison  of  Tables  10  and  12. 

- . 6  Examples 

To  illustrate  the  use  of  these  prediction  methods,  we  consider  some  examples  of 
the  sort  of  problem  to  which  they  may  be  applied. 

(i)  For  constellations  of  satellites  in  12-hour  circular  orbits,  with  a  minimum 
working  elevation  angle  at  the  Earth's  surface  of  5°  (ie  $  71.2°),  what  are 

the  numbers  of  satellites  likely  to  be  necessary  to  provide  continuous  multiple  whole- 
Earth  coverage  to  the  following  standards:  (a)  four-fold,  (b)  six-fold  and  (c)  eight-fold 
continuous  coverage? 

L'sing  values  of  c~  from  Table  15  for  n  -  4,  6  and  8,  we  obtain  the  following 
n 

values  of  T  : 

(a)  for  n  -  4,  T  =  14.18,  15,75,  18.24,  15.20,  17.58  and  21.90; 

(b)  for  n  =*  6,  T  -  20.25,  22.07,  24.93,  21.43,  24.17,  28.86;  and 

(c)  for  n  -  8,  T  -  25.90,  27.83,  30.77,  27.17,  29.98,  34.65. 

tv'e  may  interpret  these  as  meaning  that,  if  we  are  primarily  concerned  to  provide  the 

specified  coverage  with  the  minimum  number  of  satellites  and  are  not  concerned  if  the 

value  of  D„.  is  small,  we  are  likelv  to  be  able  to  meet  this  objective: 

Min 

(a)  for  four-fold  coverage,  with  between  15  and  18  satellites,  with  16  as  the 
most  likely  number  required; 

(b)  for  six-fold  coverage,  with  between  21  and  24  satellites,  with  22  as  the 
most  likely  number  required;  and 

(c)  for  eight-fold  coverage,  with  between  26  and  30  satellites,  with  28  as  the 
most  likely  number  required. 

(In  practice,  Table  10  shows  chat  15  satellites  at  ;  just  meet  the  requirement  for 

op  t 

four-fold  coverage,  and  that  24  satellites  meet  the  requirement  for  six-fold  coverage  with 
a  substantial  margin,  suggesting  that  a  smaller  number  would  be  adequate.) 

If,  on  the  other  hand,  we  are  concerned  that  while  meeting  the  coverage  requirement 
we  should  also  have  as  large  as  possible  a  value  of  »  ue  are  likely  to  be  able  to 

meet  this  objective: 

la)  for  four-fold  coverage,  wich  becween  16  and  21  satellites,  with  18  as  the 
most  likely  number  required; 

lb)  for  six-fold  coverage,  with  between  22  and  28  satellites,  with  24  as  the 
most  likely  number  required;  and 

(c)  for  eight-fold  coverage,  wich  becween  28  and  34  satellites,  with  30  as  the 
most  likely  number  required. 


(In  practice.  Tables  12  and  13  suggest  that  16  satellites  at  ^  might  just  meet  the 

requirement  for  four-fold  coverage,  and  24  satellites  that  for  six-fold  coverage.) 

(ii)  If  we  are  interested  in  eight-fold  coverage  with  a  large  value  of  3..,,.  , 

(a)  what  values  of  D  v  are  likely  to  be  obtainable  and  (b)  what  appear  to  be  the  most 
MIN 

promising  patterns  for  inicial  investigation,  preferably  avoiding  patterns  having  a  large 
value  of  P  ? 


(a)  From  Table  16,  the  predicted  largest  value  of  for  30  satellites  is 

28.1°,  the  upper  limit  associated  with  28  satellites  is  29.9°  and  the  lower 

limit  associated  with  34  satellites  is  25.3°.  Somewhat  smaller  values  would 

result  if  we  used  the  pattern  giving  the  second-largest  or  third-larges t  value,  or 

used  an  inclination  other  than  5  .  . 

opt.D 

(b)  The  most  likely  pattern  to  provide  the  required  conditions  is  the  30-satellite 
pattern  in  the  4:3  series;  from  Table  17  this  is  pattern  30/10/6.  The  next  most 
likely  patterns  are  30/15/5  and  30/15/12.  However,  considering  the  preference  for 
relatively  few  plane3  (small  value  of  P  ),  it  would  be  worth  investigating  28/7/2 
and  30/6/0;  or  if  none  of  these  proved  satisfactory,  32/8/3  might  be  suitable. 

(iii)  What  is  the  likely  value  of  Rj_^  g  for  28  satellites? 

Using  values  of  1/a  for  n  ■  8  from  Table  14,  the  likely  value  (at  :  )  is 

n  opt 

71.0°,  with  maximum  and  minimum  values  of  74.6  and  68.5  .  The  likely  value  of 

R,  „  at  5  _  is  73.7°,  with  maximum  and  minimum  values  of  79.2°  and  70.1°. 

Max, 8  opt.D 

5  USE  OF  Dmin  (OR  \  )  AS  SOLE  CRITERION  IN  A.  PRELIMINARY  STUDY  Of  NON-UNIFORM 
ORBITAL'  PATTERNS 

5 . 1  Background 

The  particular  numerical  prediction  methods  discussed  in  sections  4.3  and  a. 4  were 
derived  only  for  delta  patterns.  However,  the  conclusions  developed  in  sections  4.3  and 
4.4  regarding  the  desirability  of  identifying  those  patterns  providing  the  largest  values 
of  (°r  '  )  might  be  expected  to  have  somewhat  wider  applicability,  and  the 

simplicity  of  calculating  inter-satellite  distances  makes  this  appear  a  suitable  criterion 
for  use  in  preliminary  studies,  including  drawing  up  shorc-liscs  of  patterns  for  subse¬ 
quent  detailed  examination  against  the  particular  requirements  and  criteria  developed  for 
the  system  under  consideration.  provides  a  fully  adequace  comparison  between 

patterns  containing  similar  numbers  of  satellites,  but  1  provides  a  better  guide  to  the 
standards  being  achieved,  in  comparison  with  the  best  which  might  be  expected,  when 
different  numbers  of  satellites  are  involved. 

^MIN  wa3  necessitY  used  as  the  sole  criterion  in  the  brief  study  of  non-uniform 
patterns,  and  their  comparison  with  18-satellite  delta  patterns.  This  study  was 
occasioned  by  proposed  changes  to  plans  for  the  OPS  satellite  navigation  system;  the 
writer  did  not  have  access  to  a  computer  program  providing  values  of  Che  criceria  (GDOP 
and  PDOP)  normally  used  for  that  system,  and  while  ic  mighc  have  been  desirable  to  make 
use  of  the  full  or  shortened  version  of  the  RAE  program  COCO,  initially  developed  in  1971 
and  used  until  1976  to  calculate  the  values  of  R^  initially  published  in  Refs  3 


j-c  it  %■  .1.5  r. oc  feasible  in  the  short  tine  then  available  to  adapt  this  program  tor 
.ir.alvsis  non-uniform  patterns  'and  it  would  also,  with  such  patterns,  have  required 
sues cancial ly  n.ore  computer  tine).  However,  the  snail  portion  of  che  program  which  was 
used  to  calculate  values  of  3V...  could  be  adapted  relatively  quickly,  and  so  was  used 
in  this  study;  it  provides  a  good  example  of  what  may  be  achieved  using  Dv  as  sole 
cri terion. 

Appendix  A  of  Ref  4  discussed  the  application  of  the  delta  pattern  analysis  to 
satellite  navigation  syscems  such  as  CPS;  showed  how,  among  24-satellite  three-plane 
patterns,  the  pattern  24/3/2  appeared  preferable  to  the  pattern  24/3/1  which  seemed  from 
one  account  to  have  been  selected  ( though  it  was  subsequently  confirmed  that  24/3/2 
was  indeed  che  chosen  pattern);  and,  considering  partial  patterns  of  18  instead  of  24 
satellites,  indicated  that  18/3/0  and  18/6/2  appeared  to  be  the  best  three-plane  and 
six-plane  18-satellite  delta  paccerns  respectively,  with  18/6/2  probably  superior. 

In  1980  it  became  known  that  it  was  intended  to  reduce  the  size  of  the  initial 
operational  GPS  constellation  from  24  satellites  Co  18,  using  12-hour  orbits  of  55° 
inclination,  but  that  it  was  proposed  to  retain  the  option  to  build  up  later  co  the 
originally  planned  number  of  24  satellites.  It  was  possible  that,  in  addition  to  the  13 
active  satellites  in  orbit,  three  further  satellites  might  be  placed  in  orbit  as  inactive 
spares.  It  would  therefore  be  desirable  co  identify  patterns  which  could  meet  the 
requirements  using  only  18  satellites,  could  be  expanded  to  include  24  satellites,  and 
could  accommodate  21  satellites  in  a  suitable  manner. 

The  printouts  from  the  runs  of  the  computer  program  COCO  covering  18-satellite 
delta  patterns  which  had  been  made  in  late  1974  and  early  1975  were  re-examined  for 
further  information  relevant  co  che  choice  of  a  system  to  operate  at  55°  inclination. 

In  : 2-hour  orbits  the  value  of  ^  corresponding  Co  zero  elevation  is  76.1°, 

w:ule  cnac  corresponding  Co  5°  elevation  i.  71.2  .  For  GPS  four  well-spaced 
satellites  must  be  visible  above  5°  elevation  at  all  times,  so  in  drawing  up  a  short¬ 
list  of  potentially  suitable  patterns  it  was  considered  essential  Co  have  a  value  of 

R,  not  exceeding  71.2°,  and  the  following  arbitrarv  limits  were  also  set: 

Max,-  -  ' 

R.,  ,  not  exceeding  50  ,  R,,  .  not  exceeding  85  ,  and  D„.  not  less  than 

Max,!  6  >lax,o  6  Mm 

10'".  Seven  13-satellite  patterns  were  found  Co  meet  all  these  criteria  at  55 
inclination,  as  shown  below: 


less  Chan  76.1  , 


Ic  was  noted  that  18/13/6  and  18/18/16  both  have  values  of  8^  5 

showing  that  wich  these  pacterns  cnere  would  never  be  less  than  five  satellites  above  zero 
elevation,  while  13/6/:,  13/9/6  and  13/18/2  have  values  only  slightly  greater  chan 
76.1°,  showing  Chat  any  disappearances  of  a  fifcn  satellite  would  be  rare  and  brief. 
Patterns  18/9/6,  18/6/2  and  18/18/16  have  the  largest  values  of  D,^  •  However,  of  the 

Seven  patterns  only  18/3/0  and  18/6/2  could  be  expanded  to  26-satellite  delta  patterns; 
.3/3/0  could  also  be  expanded  to  a  21-sacellite  delta  pattern,  whereas  three  spares  added 
to  the  six-plane  pattern  18/6/2  could  be  placed  only  in  alternate  planes. 

A  fuller  comparison  of  che  possible  expansions  of  these  two  patterns,  based  on  the 
1976-5  computer  runs,  is  as  follows: 


Pattern 

Slax,  1 

*Slax,4 

^Max , 5 

^Max.b 

*Slax,7 

DMin 

18/3/0 

44.8° 

65.8° 

81.9° 

90.0° 

90.0° 

20.8° 

0 . 404 

21/3/2 

40.4° 

59.0° 

76.9° 

11.5° 

90.0° 

14.80 

0.31  1 

24/3/2 

38.  1° 

55.40 

68.2° 

69.1° 

85.2° 

10.8° 

0.242 

18/6/2 

39.70 

66.8° 

76.4° 

86.70 

90.0° 

33.3° 

0.647 

24/6/1 

38.70 

57.6° 

69.5° 

70.0° 

74.5° 

20.8° 

0.467 

This  shows  chat  che  value  of  R.  ,  for  18/6/2  is  verv  similar  to  the  values  of  R,, 
and  6  for  21/3/2;  that  boch  24/3/2  and  24/6/1  would  always  have  at  least  six 

satellites  above  5°  elevation,  but  that  24/6/1  would  always  have  at  least  seven 
satellites  above  zero  elevation,  whereas  24/3/2  would  noc;  and  chat  24/6/1  has  a 
considerably  better  value  of  D  (and  )  than  24/3/2.  Despite  the  advantage  of  che 
three-plane  series  in  terms  of  the  convenient  positioning  of  spares,  the  six-plane  series 
appears  preferable  overall. 

Fig  12  illustrates  the  satellite  distributions  in  some  of  these  patterns.  Fig  12a 
shows  pattern  24/3/2  (caking  account  of  both  cross  and  circle  symbols;  the  reason  for 
using  cwo  symbols  will  be  explained  later),  while  Fig  12b-d  shows  pacterns  13/3/0,  18/6/2 
and  24/6/1  respectively;  che  more  even  distribution  of  satellites  in  the  six-plane 
patterns  is  evident  from  chese  figures. 

Initial  approaches  to  this  problem  in  the  USA  have  been  described  by  3ook  si'*”* 
and  by  Jorgensen'^.  The  obvious  first  approach  was  to  retain  a  chree-plane  configuration 
but  with  che  optimum  uniform  distribution  of  six  instead  of  eight  satellites  in  each 
plane,  ie  to  change  from  delta  pattern  24/3/2  to  delta  pattern  13/3/0.  However,  though 
this  gave  fairly  uniform  coverage,  the  performance  provided  in  terms  of  che  navigation 
criterion  PDOP  was  poor. 

Book  ez  si  then  investigated  other  distributions  of  18  satellites  in  three  planes, 
and  identified  a  non-uniform  pattern  which  provided  better  performance  in  terms  of  the 
navigation  criterion  PDOP;  this  corresponded  to  the  original  delta  pattern  24' 3/ 2  with 
two  satellites  removed  from  each  of  the  three  planes  and  the  remainder  left  in  their 
original  positions.  A3  noted  in  Appendix  A  and  Table  7  of  Ref  4,  in  12-hour  orbits 
pattern  24/3/2  would  follow  eight  separate  Earth-tracks ,  while  18/3/0  would  follow  six 


and  18/6/2  would  follow  nine;  the  six  satellites  removed  from  24/3/2  would  correspond  to 
those  following  two  adjacent  tracks,  leaving  six  tracks  only.  This  non-uniform  18- 
satellite  three-plane  pattern  is  illustrated  by  the  circle  symbols  in  Fig  12a,  with  the 
crosses  representing  the  satellites  removed  from  pattern  24/3/2.  It  may  be  noted  that 
the  values  of  D^  at  55°  for  chese  patterns  are;  24/3/2  -  10.8°  (X  •  0.242); 

18/3/0  -  20.8°  (X  *  0.404);  while  the  non-uniform  18-satellite  three-plane  pattern  has 
the  same  value  of  as  24/3/2  from  which  it  is  derived,  ie  10.8°  (X  -  0.210). 

The  orders  of  merit  given  by  PDOP  and  X  are  therefore  contradictory  in  this  case;  this 
may  be  partly  explained  by  the  fact  that  X  deals  with  worst-case  conditions,  whereas 
PDOP  is  statistical  in  nature  and  can  overlook  isolated  adverse  cases.  When  a  pattern 
is  uniformly  'good',  it  is  likely  to  be  recognised  as  such  by  both  X  and  PDOP;  whereas, 
given  a  choice  becween  one  which  is  uniformly  'mediocre'  and  one  which  is  sometimes  'good' 
and  sometimes  'poor',  the  former  will  be  favoured  by  X  and  the  latter  by  PDOP. 

Another  non-uniform  18-satellite  constellation  was  subsequently  proposed  by 
13  14 

Blake  '  .  In  Ref  4  it  was  remarked  that  other  authors  had  apparently  not  recognised 

single-satellite-per-plane  patterns  as  important  members  of  the  family  of  delta  patterns; 

I  o 

these  were  subsequently  discussed  by  Ballard  under  the  title  of  'rosette  constellations’, 

and  Blake's  proposal  was  for  a  'modified  rosette’  in  which  six  evenly-spaced  satellites 

were  removed  from  the  delta  (or  rosette)  pattern  24/24/2  (ie  one  from  every  fourth 

plane)  and  the  phase  spacing  (but  not  the  longitude  spacing)  of  the  remainder  was  then 

adjusted  to  be  the  same  as  that  for  pattern  18/18/2,  ie  phase  difference  between 

adjacent  planes  changed  from  30°  to  40°.  Pattern  24/24/2  (with  the  six  satellites 

removed  indicated  by  crosses  and  those  remaining  indicated  by  circles)  and  Blake's  18- 

satellite  modified  rosette  derived  from  it  are  illustrated  in  Fig  I3a4b  respectively. 

Jorgensen1’’  has  indicated  that,  in  terms  of  percentage  visibility  and  cumulative 

accuracy  probability,  the  modified  rosette  gives  slightly  better  results  than  the  non- 

uniform  three-plane  pattern,  but  that  pattern  18/6/2  appears  superior  to  both.  It  may 

be  noted  chat  D„.  for  pattern  24/24/2,  and  for  the  18-satellite  pattern  left  when  six 
Min  Q 

satellites  are  removed  from  it,  is  12.0  (X  »  0.233  for  18  satellites);  for  Blake's 
modified  rosette  with  adjusted  spacing  it  is  18.4°  (X  »  0.358);  and  for  pattern 
18/6/2  it  is  33.3°  (X  «  0.647).  The  X  values  therefore  agree  with  Jorgensen's  order 
of  merit  in  these  cases. 

It  appeared  to  the  present  writer  that  31ake's  proposal  might  be  capable  of  further 
development;  this  led  to  the  brief  study7  whose  results  are  described  in  the  following 
section. 

5 . 2  Non-uniform  patterns 

Some  larger  delta  patterns,  involving  a  total  number  of  satellites  vT)  which  is  not 
a  prim  number,  may  be  regarded  as  being  built  up  from  several  similar  delta  sub- 
patterns,  evenly  spaced;  and  some  patterns  which  are  not  themselves  delta  patterns  may 
nevertheless  also  be  regarded  as  being  built  up  from  several  similar  de lea  sub-patterns, 
chough  with  chese  not  all  evenly  spaced.  Thus  pattern  2-,  24/2  may  be  regarded  as  built 
up  of  four  examples  of  patcern  6/8/2,  spaced  by  15°  in  longitude  and  by  30°  in 
phase.  Blake's  modified  rosette  consists  of  three  examples  of  pattern  6/6/2,  spaced  by 
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15°  in  longitude  (with  one  30°  gap)  and  by  40°  in  phase;  clearly  in  such  a  case 
the  phase  spacing  of  the  sub-patterns  could  cake  any  suitable  value.  In  the  initial 
account^  of  this  study  che  name  'omega  patterns'  was  adopted  for  such  non-uniform 
patterns  built  up  from  delta  sub-patterns;  and,  considering  specifically  those  which  have 
evenly-spaced  omissions  from  an  otherwise  even  spacing  of  planes  in  a  rosette  pattern, 
an  omega  pattern  identification  code  T/P/F/W  was  adopted,  based  on  the  delta  pattern 
identification  code  T/P/F,  which  is  retained  for  the  delta  sub-pattern  in  the  form 
Tg/Pg/F,_.  In  the  omega  pattern  code  T/P/F/W,  T  is  the  actual  number  of  satellites  in 
the  pattern  (after  any  deletions);  P  is  che  total  number  of  planes,  both  occupied  and 
empty,  corresponding  to  the  even  spacing  of  the  occupied  planes;  F  is  the  same  as  the 
value  of  F^  for  the  sub-pattern;  and  W  ,  expressed  in  units  of  360°/P,  is  the  phase 
spacing  between  the  sub-patterns,  measured  in  che  same  sense  as  F  ,  and  can  cake  any 
real  value  from  0  to  P  (with  W  «  P  equivalent  to  W  »  0).  Note  that  P^  »  T<.  is  equal 
to  the  value  of  (P  -  T)  for  the  omega  pattern,  which  must  be  a  factor  of  both  T  and  P  . 

Using  this  identification  code,  removing  one  of  the  four  6/6/2  sub-patterns  from 
che  delta  pattern  24/24/2  leaves  the  omega  pattern  18/24/2/2.0,  while  changing  the  phase 
spacing  between  che  sub-patterns  from  30°  to  40°  in  accordance  with  Blake's 
proposal  produces  che  omega  pattern  18/24/2/2.667. 

Clearly  the  characteristics  of  an  omega  pattern  are  partly  determined  by 'those  of 

the  delta  sub-pattern  from  which  it  is  built  up  and,  in  particular,  the  value  of  Dv,^n 

for  an  omega  pattern  cannot  exceed  that  for  the  sub-pattern  at  the  same  inclination. 

Hence,  if  the  sub-pattern  has  D„.  *  0°,  any  omega  pattern  based  on  it  will  also  have 

run 

»  0°;  and  it  appears  unwise  to  choose  a  sub-pattern  having  a  value  of  Dj^n 

substantially  less  than  that  of  the  delta  pattern  giving  the  largest  value  of  D,,.  for 

Min 

the  value  of  T  to  be  used  for  the  omega  pattern  ( eg  38.1°  for  18  satellites),  since 
the  sub-pattern  would  chen  probably  be  the  limiting  factor  on  the  value  of  Dj^  for 
the  omega  pattern. 

Considering  all  possible  24-satellite  rosette  patterns  composed  of  6-satellite 
rosette  sub-patterns,  we  find  that  there  are  six  possible  sub-patterns  which  between 
them  produce  24  rosette  patterns,  each  of  which  leaves  an  omega  pattern  when  one  sub¬ 
pattern  is  deleted;  these  are  identified  as  follows: 


Sub-pattern 

VVFs 

Rosette  pattern 
T/P/F 

Omega  pattern 
T/P/F/W 

Sub-pattern 
T  /P  /F 

s'  s  s 

Rosette  pattern 
T/P/F 

Omega  pattern 
T/P/F/W 

6/6/0 

' 

”  24/24/0 
24/24/6 
24/24/12 
24/24/18 

18/24/0/0.0 

18/24/0/6.0 

18/24/0/12.0 

18/24/0/18.0 

6/6/3 

f  24/24/3 

J  24/24/9 
]  24/24/15 
l  24/24/21 

18/24/3/3.0 

18/24/3/9.0 

18/24/3/15.0 

18/24/3/21.0 

6/6/1 

1 

f  24/24/1 
24/24/7 
24/24/13 
,  24/24/19 

18/24/1/1.0 

18/24/1/7.0 

18/24/1/13.0 

18/24/1/19.0 

6/6/4 

f  24/24/4 

J  24/24/10 
]  24/24/16 
l  24/24/22 

18/24/4/4.0 

18/24/4/10.0 

18/24/4/16.0 

18/24/4/22.0 

6/6/2 

'  24/24/2 
24/24/8 
24/24/14 
24/24/20 

18/24/2/2.0 

18/24/2/8.0 

18/24/2/14.0 

18/24/2/20.0 

6/6/5 

(  24/24/5 

J  24/24/11 
]  24/24/17 
l  24/24/23 

18/24/5/5.0 

18/24/5/11.0 

18/24/5/17.0 

18/24/5/23.0 

~h^s%  for  che  crr.ega  patterns,  when  ve  come  to  consider  values  of  0^^  over  Che  full 

ranee  of  values  of  T*  from  0  to  -3.999,  ve  know  chat  for  each  value  of  F  there  are 

four  particular  values  of  W  at  which  che  value  of  DXi  will  be  the  same  as  for  the 

. .  l  n 

corresponding  rosette  pattern  and  that,  for  all  values  of  V  ,  D„.  will  not  exceed  the 
r  “  Mm 

value  for  the  sub-pattern. 

From  section  3.5  of  Ref  4  we  know  that  a  delta  pattern  will  have  D,,.  *  0°  if  boch 

r  Mm 

?  and  (T/P  -  F)  are  even  numbers.  Here  this  applies  to  the  delta  sub-patterns  6/6/1, 
6/6/3  and  6/6/5,  and  to  all  che  rosetce  and  omega  patterns  based  upon  them;  these  may 
therefore  be  excluded  from  further  discussion. 

For  sub-pattern  6/6/2  we  have  already  (in  the  previous  section)  noted  values  of 

D,,.  (at  35°  inclination)  at  W  »  2.0  and  2.667;  the  values  at  W  »  3.0,  14.0  and  20.0 
Mm 

correspond  to  those  for  the  ocher  three  rosette  patterns,  and  we  may  calculate  further 

values  of  D,,.  at  intermediate  values  of  W  to  obtain  the  complete  picture  of  the 
Mm  oo 

variation  of  Dw.  with  W  .  We  find  chat  falls  from  6.1  at  W  *  0  to  4.4  at 

Mm  Mm 

1.0,  then  rises  to  12.0°  ac  2.0,  15.2°  at  2.333,  18.4°  at  2.667,  and  21.3°  (the  limiting 
value  for  sub-pattern  6/6/2)  at  3.0.  It  then  falls  to  18.2°  at  W  »  3.333,  6.1°  ac  4.0, 
12.2°  ac  5.0,  0.3''  ac  6.0,  and  rises  again  to  the  limiting  value  of  21.3°  at  7.0,  7.25 
and  7.3,  before  falling  to  19.7°  ac  3.0,  11.2°  at  9.0  and  0.3°  at  10.0.  Similar  fluctua¬ 
tions  of  D,,.  continue  across  the  rest  of  the  range  of  values  of  W  .  Thus  we  see  chat 
Mm 

we  can  improve  on  the  value  of  D^; ^  of  18.4  ac  the  value  of  K  proposed  by  Blake, 
reaching  che  limiting  value  of  21.3°  for  values  of  K  in  the  neighbourhood  of  3.0  and 
7.25.  Values  of  POOP  would  not  necessarily  vary  in  similar  manner,  buc  their  behaviour 
would  appear  worth  investigating.  Satellite  positions  in  chese  patterns  are  lisced  in 
Table  18,  and  pattern  18/24/2/7.25  is  shown  in  Fig  13c. 

Sub-patterns  6/6/0  and  6/6/4,  on  che  other  hand,  have  values  of  D^n  of  33.3°  and 

70.0°  respectively  ac  55°  inclination;  there  is  therefore  some  prospect  that  omega 

patterns  based  on  either  of  chese  might  prove  superior  to  those  based  on  sub-pattern 

6/6/2.  A  limited  investigation  was  made  of  omega  patterns  based  on  sub-pattern  6/6/4. 

From  the  rosetce  patterns,  we  find  that  D,^  *  6.1°  ac  V  *  4.0,  2.4°  at  10.0,  6.1°  at 

16.0  and  21.2°  ac  22.0;  near  the  latter  value  of  K  ,  it  falls  to  18.3°  at  22.2,  but 

rises  to  24.2°  at  21.8  and  Co  a  peak  of  28.1°  at  W  »  21.541  before  falling  to  24.0°  at 

21. a,  12.2°  at  21.0  and  2.4°  at  20.5,  then  rising  to  a  smaller  peak  at  26.8°  ac 

W  »  19.667.  Thus  a  0.,.  of  at  least  28.1°  can  be  obtained  bv  changing  to  a  different 
Mm 

sub-pattern.  Satellite  positions  in  these  patterns  are  lisced  in  Table  19,  and  pattern 
I  S/24/4/2 1  . 341  is  shown  in  Fig  13d.  This  is  particularly  interesting  in  that  it  shows 
that,  in  optimising  for  D,,.  ,  we  have  arrived  ac  a  2-*-plane  pattern  (shown  bv  circles) 

which  matches  very  closely  (though  not  exactly)  the  delta  pattern  18/18/6  (shown  by  *s) 
which  was  one  of  the  short-list  of  seven  whose  charac ter i s tics  were  quoted  in  the 
previous  section. 

We  have  thus  far  followed  31ake  in  assuming  that  we  should  look  for  an  13-satellite 
omega  pattern  expandable  to  a  2--sacelliCe  rosetce  pattern.  However,  there  are  some 
advantages,  in  the  circumstances  described  in  che  preceding  section,  in  considering  an 
'.6-satellite  omega  pattern  based  on  a  21-satellite  rosette  pattern,  is  in  optimising  the 
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system  at  the  21-satellite  level  rather  than  at  the  24-satellite  level.  A  2l-satellice 
delta  pattern  may  be  considered  to  be  composed  of  seven  3-satellite  sub-patterns;  to  meet 
the  requirement  under  discussion,  any  six  of  these  sub-patterns  may  be  active  and  any  one 
the  inactive  spare  set.  Different  sub-patterns  could  be  designated  as  the  spare  set  in 
turn  at  di£fe..ent  times,  so  equalising  satellite  life  expectancy  over  che  whole  system, 
and  when  any  one  satellite  failed  the  set  of  uhich  it  formed  part  could  be  designated  as 
the  spare  set  until  a  replacement  had  been  launched,  thus  completely  'hiding'  a  single 
satellite  failure  (though  such  changes  would  require  some  re-adjustment  of  phasing  among 
the  accive  satellites  if  the  value  of  W  did  noc  correspond  directly  to  a  rosette 
pattern).  Hence,  though  a  21-plane  18-satellite  pattern  provides  a  less  satisfactory 
basis  than  a  24-plane  18-sacellite  pattern  for  possible  later  expansion  to  a  24-satellite 
system,  since  this  would  have  co  be  non-uniform,  the  2l-plane  constellation  seems  likely 
to  provide  a  more  failure-resistant  system  at  che  more  critical  initial  stage  when  only 
IS  satellites  would  be  active,  with  three  inactive  spares. 

Of  che  21-plane  delta  (rosette)  patterns  which  might  form  a  basis  for  such  a 
system,  che  two  most  promising  are  probably  21/21/9  and  21/21/19.  Ac  55°  inclination 
21/21/19,  and  hence  also  the  omega  pattern  18/21/1/19.0,  has  »  24.2°;  this  is 

increased  as  W  is  reduced,  wich  a  maximum  of  D .  »  28.1°  for  pattern  18/21/1/18.772. 

Ac  che  same  inclination  21/21/9,  and  hence  also  che  omega  pattern  18/21/0/9.0  (see 
Table  20  and  Fig  13e),  has  »  30.4°.  This  is  at  least  a  local  maximum,  falling 

as  W  is  either  increased  or  decreased  from  the  value  W  »  9.0;  the  satellites  would  not 
(since  U  is  unchanged  from  the  rosette  pattern)  need  phase  adjustments  when  sub-patterns 
changed  between  active  and  inactive  status;  and  chis  local  maximum  may  be  the  largest 
value  of  obtainable  wich  a  21-plane  18-satellite  pattern  (chough  noc  all  values  of 

W  over  the  full  range  from  0  co  20.999  have  been  checked  to  confirm  this).  Indeed,  the 
example  of  pattern  18/24/4/21.541  raises  che  possibility  that  a  21-plane  pattern  might  be 
found  which  would  closely  match  the  delta  pattern  18/9/6,  though  this  has  not  been 
pursued  further. 

Thus  the  largest  values  of  D„.  for  18-satellite  constellations  at  55°  inclina- 

Mtn 

cion  which  we  have  found  during  this  brief  study  using  the  non-uniform  omega  patterns  are 
30.4°  for  che  21-plane  pattern  18/21/0/9.0  and  28.1°  for  the  24-plane  pattern 
18/24/4/21.541.  These  values  compare  with  values  for  che  uniform  delca  patterns  of  33.3° 
for  che  six-plane  pattern  18/6/2,  28.2°  for  the  18-plane  pattern  18/18/16,  and  20.8°  for 
Che  three-plane  pattern  18/3/0.  Of  these  pacterns,  positioning  of  three  inactive  spares 
and  bringing  chem  into  service  after  a  failure  appears  to  be  mosc  readily  accomplished  in 
patterns  18/21/0/9.0,  18/6/2  and  18/3/0,  while  expansion  co  24  satellites  would  be  mosc 
readily  performed  wich  patterns  18/24/4/21.541,  18/6/2  and  18/3/0.  Overall,  delta  paccern 
18/6/2  may  well  be  considered  to  have  mosc  advantages  (though  some  other  relevant  factors, 
such  as  launching  arrangements,  have  not  been  considered  here);  Brady  and  Jorgensen1 
indicated  that,  in  1981,  this  pattern  was  (subject  to  some  further  consideration  of 
possible  advantages  of  a  three-plane  pattern  as  regards  Shuttle  launching)  che  current 
baseline  orbital  configuration  for  the  operational  Phase  III  GPS. 
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This  Report  has  described  several  separate,  though  related,  studies  performed 
during  the  period  1 977-1 9S2.  The  principal  conclusions  emerging  from  these  studies  are 
as  follows: 

ill  An  examination  of  Che  coverage  characteristics  of  regular  tesselations  of  an 
infinite  plane  and  of  regular  polyhedra  has  led  to: 

la)  the  identification  of  trends,  with  respect  to  numbers  of  satellites  and 
degrees  of  coverage,  which  appear  applicable  also  to  practical  constellations  of 
satellites  in  Earth  orbit; 

(b)  the  use  of  alternative  versions  (l/o  and  '•),  normalised  bv  reference  to  the 

n 

characteristics  of  an  ideal  triangular  tesselation,  of  the  criteria  (R^„  ^  and 
D  )  which  have  been  used  in  previous  coverage  studies;  and  hence 

i.c)  che  fornulacion  of  methods  of  predicting  approximate  values  of  these  criteria 
which  may  be  expected  to  be  achievable  with  larger  numbers  of  satellites  and  higher 
degrees  of  coverage  than  have  yet  been  examined  in  detail  for  practical  satellite 
conste  nations . 

■  3)  An  attempt  to  produce  a  combined  criterion  incorporating  both  and  D  , 

has  demonstrated  the  much  stronger  influence  of  D  ,  and  suggests  that  it  would  be 
appropriate  to  use  •  “hich  is  much  the  simpler  to  calculate,  as  the  sole  criterion 

in  drawing  up  a  short-list  of  orbital  patterns  for  detailed  study  against  any  particular 
coverage  requirement. 

(3)  It  has  been  confirmed  that  the  delca  patterns  giving  the  largest  values  of 
are  those  which  produce  single  non-self-crossing  Earth-tracks  at  certain  altitudes; 
these  patterns  can  be  identified  by  methods  which  were  developed  in  a  previous  report-*. 

i-)  An  examination  of  non-uniform  patterns,  using  as  criterion,  has  suggested: 

ia)  chat  ic  is  possible  to  improve  on  results  obtained  elsewhere  by  choosing  che 
most  suitable  sub-pattern  and  then  optimising  the  phasing  between  adjacent  planes, 
and 

lb)  chat  in  a  non-uniform  syscem  including  both  active  satellites  and  spares  in 
orbit,  ic  may  be  desirable  to  optimise  for  che  total  number  of  satellites  in  orbit. 

ij)  Results  obtained  in  section  5  have  been  discussed  with  reference  Co  che  require¬ 
ments  of  the  GPS  navigation  satellite  syscem.  Subject  to  satisfactory  launching  arrange¬ 
ments  being  determined,  it  appears  chat  the  delta  pattern  13/6/2,  to  which  attention  was 
originally  drawn  in  Appendix  A  of  Re:  A,  may  well  have  most  overall  advantages  among 
potential  lo-satellice  cons te  11  at  ions  ,  and  ic  is  understood  to  have  been  selected  by  the 
programme  authorities  as  the  CPS  baseline. 
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Table  5 


DERIVATION  OF  I/;  :  REGULAR  TESSELATION'S 

n 


RXAX,n’”L 

“n 

1  tz 

n 

■ 

s 

H 

T 

S 

H 

T 

s 

H 

I 

0.57'A 

0.7071 

1.0000 

1.0997 

1.2533 

1.5551 

0.9093 

0.7979 

0.6-30 

2 

I .3000 

1 . 0000 

1 .0000 

1.3468 

1.2533 

1.0996 

0.7-25 

0.7979 

0.9094 

3 

1  .0000 

1.1180 

1 .3229 

1 .0996 

1.14-1 

1.1878 

0.9094 

0.8741 

0.8-19 

4 

1.13-7 

1.1785 

1 .-000 

1 .0996 

1 • 0444 

1 .0886 

0.9094 

0.9597 

0.9136 

5 

1.3::? 

1 .581  1 

1.7321 

1.1268 

1.2533 

1 . 2046 

0.8875 

0.7979 

0.8301 

3 

i . 3-7; 

1.5811 

1 .7321 

1 .0475 

1.1441 

1.0997 

0.9546 

0.8741 

0.909- 

7 

1.5276 

1.5311 

2.0000 

1.0997 

1.1756 

0.9093 

0.9441 

0.8507 

3 

1.7321 

1 .6667 

2.0000 

1.1664 

1.0444 

1.0996 

0.8574 

0.957- 

0.9094 

9 

1.7321 

1 .8028 

2.0817 

1  .0997 

1 .0651 

1 .0791 

0.9094 

0.9389 

0.9267 

10 

1.7500 

2.0000 

2.2432 

1 .0540 

1 .1210 

1.1041 

0.9488 

0.892) 

0.9057 

1 1 

1.8359 

2.0156 

2.3065 

1 .0543 

1.0772 

1.0815 

0.9485 

0.9284 

0.9247 

1  2 

1.  8571 

2. 102- 

2.3333 

1.021! 

1.0757 

1.0475 

0.9794 

0.9296 

0.9547 

13 

2.0817 

2. 1 250 

2.6458 

1 .0997 

1 .0446 

1.1412 

0.9094 

0.9573 

0.8763 

14 

2.0817 

2.2361 

2.6458 

1 .0597 

1.0593 

1.0997 

0.9437 

0.94-1 

0.909- 

1  5 

2.1794 

2.3049 

2.6458 

1.0718 

1 .0548 

1 .0624 

0.9330 

0.9480 

0.9413 

IS 

2. 132- 

2.3570 

2.6458 

1 .0392 

1 .0444 

1 .0286 

0.9623 

0.9575 

0.9722 

i  7 

2.2-52 

1 .037! 

0.9642 

IS 

2.3065 

1 .0354 

0.9658 

i  9 

2.3333 

1.0195 

0.9808 

20 

2.6-53 

1.1268 

0.3875 

2  ! 

2 .  Su»38 

1 .0997 

0.9094 

2  2 

2 .6-58 

i  .  0  7 44 

0.9308 

23 

2 . 6458 

1 .0508 

0.9517 

2- 

2.6-58 

1 .0286 

0.9722 

2  5 

2 . 6674 

1 .0161 

0.9342 

r 

1 .0000 

1 .0000 

I .0000 

: .0000 

1 .0000 

1 .0000 

T  »  triangular  cesselacion 
S  *  square  tesselation 
H  ■  hexagonal  tesselacion 
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Table  8 


DERIVATION'  OF  VALUES  OF 


^HAX,n 


FOR  REGULAR  POLYREDRA 


Polyhedron 

Radius 

(deg) 

■ 

Relevant  values  of  n 

Tetrahedron 

70.5288 

1 

109.4712 

- 

m 

2 

Octahedron 

54.7356 

3,  o 

1 

90.0000 

2,  3 

125.2644 

3,  3 

4 

Cube 

54.7356 

4.  0 

1  (2) 

70.5288 

3.  ! 

2 

90.0000 

4.  2 

3,  4 

109.4712 

3,  4 

(5) 

125.2644 

4.  4 

5,  6 

Icosahedron 

37.3774 

3.  0 

1 

63.4350 

5.  1 

2,  3  (4) 

79. 1877 

3,  3 

4 

90.0000 

4,  4 

5.  6 

100.8123 

3.  6 

(7) 

116.5650 

5.  6 

7,  8,  9 

142.6226 

3,  9 

10 

Dodecahedron 

37.3774 

5,  0 

1  (2,  3) 

41.8103 

3,  I 

*> 

54.7356 

4.  2 

3  (4) 

56.7949 

3.  3 

4 

60.7941 

3,  4 

(5) 

70.5288 

6,  4 

5  (6,  7,  8) 

79.1877 

5,  5 

6,  7,  8 

90.0000 

4,  8 

9,  10 

100.8123 

5.  10 

(11,  12,  13) 

109.4712 

6,  10 

11,  12,  13  (14) 

1 19.2059 

3,  13 

14 

123.2051 

3.  14 

05) 

125.2644 

4,  14 

15  (16) 

138.1897 

3,  16 

(17) 

142.6226 

5,  15 

16,  17,  18 

Table  9 


DERIVATION  CF  1,;  :  REGULAR  ?0 LVHEDEA 

n 


Polyhedron 

n 

^lAX.n 

(deg) 

~  n 

l/c 

n 

Tecrahedron 

■ 

70.5238 

1.2310 

nm 

H 

109.4712 

1.3510 

§§i§| 

Octahedron 

1 

54. 7356 

1.1700 

0.85-7 

6 

90. 0000 

1.3603 

0.7351 

3 

90.0000 

1,1107 

0.9003 

4 

125.2644 

1.3388 

0.7469 

Cube 

i 

54.7356 

1.3510 

0.7402 

2 

70.5288 

1.2310 

0.6124 

3 

90.0000 

1.2325 

0.7797 

4 

90.0000 

1.1107 

0.9003 

5 

125.2644 

1.3827 

0.7232 

6 

125.2644 

im 

0.7922 

Icosahedron 

1 

37,3774 

1.1299 

0.8850 

2 

63.4350 

1 .3560 

0.7375 

3 

63.4350 

1.1071 

0.9032 

4 

79.1877 

1 . 1969 

0.8355 

5 

90.0000 

1.2167 

0.3219 

6 

90.0000 

1.1107 

0.9C03 

7 

1 16.5650 

1 .3317 

0.7503 

8 

1 16.5650 

1 . 2*58 

0  •  3  C  "*  7 

9 

116. 5650 

1.1746 

0.3514 

10 

142.6226 

1.363* 

2. '335 

Dodecahedron 

1 

37.3774 

1  . 4357 

•7.6555 

41 .8103 

1.153S 

? . 566" 

3 

54. 7356 

1.2333 

v  .  w  ■ 

- 

56 . 7949 

1 .1033 

0.9023 

5 

70.5288 

1.2310 

C.312* 

6 

79.1877 

1.2617 

0.7*26 

7 

79.1877 

1 .  16S1 

0.556  1 

8 

79.1377 

1  .09  26 

0.9152 

9 

90.0CC0 

1  .  I7CS 

0 . 5  5  -  : 

10 

90 .0000 

1 . 1 107 

7.6003 

!  1 

109.4712 

1  .-65  1 

1  2 

109. -'12 

!  0  J  • 

13 

109. -.7  12 

1.13-9 

'. S-3- 

1  - 

1 19.2054 

1 

I  5 

125.26— 

1 ..-22 

-  7  - 

16 

u:.6::6 

1.3415 

: . 'i Sn 

i ; 

142.622:) 

1 .35  '  - 

IS 

1*2.6226 

1  .3,  19 

*• .  ■  n  5  - 

Table  10 


SMALLEST  VALUES  OF  R^^  n  FOR  DELTA  PATTERNS 


a 

1 

Pattern 

RMAX,n 

(deg) 

6 

opt 

(deg) 

°Min 

(deg) 

pn 

,/cu 

\ 

z 

n 

i 

5 

5/5/1 

69.2 

43.7 

60.9 

1.350 

0.741 

0.624 

0.462 

6 

6/6/4 

66 . 4 

53.1 

73.7 

1.419 

0 .705 

0.827 

0.583 

7 

7/7/5 

60.3 

55.7 

57.0 

1.392 

0.718 

0.691 

0.496 

8 

8/8/6 

56. 5 

61.9 

56.3 

1  .395 

0.717 

0.730 

0.523 

9 

9/9/7 

54.8 

70.5 

43.1 

1  .435 

0.697 

0.592 

0.413 

10 

10/5/2 

52.2 

57.1 

46.6 

1.441 

0.694 

0.675 

0.469 

1  1 

11/11/4 

47.6 

53.8 

49.0 

1  .378 

0.726 

0.745 

0.540 

12 

12/3/1 

47.9 

50.7 

26.0 

1.448 

0.691 

0.413 

0.285 

13 

13/13/5 

43.8 

58.4 

45.9 

1.378 

0.726 

0.758 

0.550 

14 

14/7/4 

42.0 

54.0 

42.5 

1.371 

0.729 

0.729 

0.531 

15 

15/3/1 

42.1 

53.5 

20.2 

1.423 

0.703 

0.358 

0.252 

2 

7 

7/7/2 

76.0 

61.8 

37.  1 

1  .241 

0.806 

0.450 

0.362 

8 

8/8/2 

71.0 

57.1 

9.6 

1.239 

0.807 

0.124 

0.100 

9 

9/3/2 

66.2 

62.1 

24.2 

1.225 

O.b  i  o 

0.333 

0.271 

10 

10/10/2 

64.  1 

61.6 

21.2 

1.251 

0.799 

0.307 

0.246 

1  1 

11/11/9 

62.0 

52.7 

44.  1 

1.269 

0.788 

0.670 

0.528 

12 

12/3/1 

56.6 

57.0 

18.2 

1.210 

0.827 

0.289 

0.239 

13 

13/13/3 

54.7 

52.8 

38.0 

1.217 

0.822 

0.628 

0.516 

14 

14/14/10 

52.4 

53.8 

16.7 

1.210 

0.827 

0.286 

0.237 

15 

15/3/1 

51.3 

55.3 

21  .5 

1.226 

0.816 

0.369 

0.301 

10 

10/10/8 

80.3 

60.0 

51.5 

1.279 

0.782 

0.746 

0.583 

1  1 

11/11/3 

74.6 

59.8 

3.9 

1.247 

0.802 

0.059 

0.048 

12 

12/4/2 

70.9 

60.0 

5.4 

1.237 

0.808 

0.086 

0.069 

13 

13/13/4 

68.0 

50.0 

27.1 

1.235 

0.810 

0.448 

0.362 

14 

14/14/4 

66.1 

47.6 

4.6 

1.246 

0.803 

0.079 

0.063 

15 

15/15/6 

63.2 

57.0 

41  .9 

1  .233 

0.811 

0.744 

0.603 

D 

13 

13/13/2 

77.1 

45.7 

29.3 

1.213 

0.824 

0.484 

0.399 

14 

14/14/4 

75.8 

69.6 

4.3 

1  .238 

0.808 

0.074 

0.060 

15 

15/15/2 

70.9 

55.7 

21  .8 

1.198 

0.835 

0.387 

0.323 

18 

18/3/0 

65.1 

57.5 

17.6 

1.205 

0.830 

0.342 

0.284 

21 

21/3/2 

58.8 

54.4 

15.5 

1.176 

0.851 

0.325 

0.277 

5 

18 

18/6/2 

73.6 

64.6 

20.5 

1.219 

0.821 

0.398 

0.327 

6 

24 

24/ 4/3 

68.3 

56 . 7 

14.1 

1.192 

0.839 

0.316 

0.265 

7 

24 

24/8/4 

75.8 

5°. 9 

24.2 

1.225 

0.816 

0.543 

0.443 

24/8/ 


75.8 


25 
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Table  12 


LARGEST  VALUES  OF  FOR  DELTA  PATTERNS 


1 

Pattern 

°MIN 

(deg) 

■ 

-op t  ,D 
(deg) 

n 

RMax,n 

(deg) 

“  n 

\/z 

n 

€ 

n 

5 

5/5/3 

82.2 

0.842 

51.8 

n 

75.5 

1 .474 

0.679 

0.571 

6 

6/6/4 

73.7 

0.828 

53.1 

n 

66.4 

1  .420 

0.704 

0.583 

7 

7/7/5 

68.3 

0.828 

58.6 

n 

60.7 

1  .401 

0.714 

0.591 

H 

81.5 

1  .331 

0.752 

0.622 

8 

8/8/6 

61.9 

0.802 

59.  1 

i 

56.9 

1 .405 

0.712 

0.571 

2 

74.0 

1.292 

0.774 

0.621 

9 

9/9/7 

57.9 

0.795 

62.7 

i 

55.1 

1  .441 

0.694 

0.552 

2 

69.7 

1  .290 

0.775 

0.617 

3 

86.  1 

1.302 

0.768 

0.611 

10 

10/5/2 

53.1 

0.770 

53.1 

1 

53.1 

1  .467 

0.682 

0.525 

2 

65.2 

1  .273 

0.786 

0.605 

1  1 

11/11/4 

50.1 

0.761 

53.1 

1 

47.7 

1.380 

0.725 

0.551 

2 

72.1 

1.476 

0.678 

0.516 

3 

78.0 

1.304 

0.767 

0.584 

4 

86.5 

1.252 

0.799 

0.608 

12 

12/6/3 

48.2 

0.764 

54.7 

i 

54.9 

1  .658 

0.603 

0.461 

2 

61.5 

1.314 

0.761 

0.582 

3 

72.6 

1.267 

0.789 

0.603 

13 

13/13/5 

47.3 

0.781 

57.6 

1 

44.  1 

1.389 

0.720 

0.562 

2 

60.4 

1.344 

0 . 7  44 

0.581 

3 

75.4 

1.370 

0.730 

0.570 

Qi 

80.6 

1.268 

0.789 

0.616 

1  u 

14/7/4 

44.5 

0.762 

57.3 

42.4 

1  .385 

0.722 

0.550 

60.0 

1.385 

0.722 

0.550 

67 . 5 

1.273 

0.786 

0.599 

77.  1 

1.259 

0.794 

0.605 

15 

15/15/6 

42.6 

0.756 

58. 1 

i 

46 . 6 

1.575 

0.635 

0.480 

2 

55.0 

1.313 

0.761 

0.576 

3 

63.4 

1  .237 

0.808 

0.611 

4 

78.5 

1.327 

0.753 

0.570 

16 

16/8/5 

41.4 

0.759 

59.9 

1 

40.8 

1.423 

0.703 

0.534 

i 

53.5 

1.320 

0.758 

0.575 

3 

65 ,  l 

1.312 

0.762 

0.579 

4 

70.8 

1.235 

0.810 

0.615 

5 

30.4 

1.25- 

0.797 

0.605 

(Cl 


Table  13 


SECOND-LARGEST  VALUES  OF  DWTV  FOR  DELTA  PATTERNS 

MIN 


Table  I 4 


VALUES  OF  \/z  FOR  USE  IN  PREDICTIONS 
n 


For  minimum 

SlAX.n 

For 

maxi mum 

dmin 

0.20 

0.26 

0.34 

0.24 

m 

0.43 

l 

0.8154 

0.7600 

0.6861 

0.7785 

Q.7046 

Q . 603 1 

2 

0.8296 

0.7784 

0.7103 

0.7955 

0.7273 

0.6336 

3 

0.8427 

0.7955 

0.7325 

0.81  12 

0.7483 

0.6617 

•4 

0.8548 

0.8112 

0.753! 

0.8257 

0.7676 

0.6878 

5 

0.8659 

0.8257 

0.7721 

0.8391 

0.7855 

0.7113 

6 

0.8762 

0.8391 

0.7896 

0.3515 

0.8020 

0.7339 

7 

0.8858 

0.8515 

0.8058 

0.8629 

0.8172 

0 . 7544 

8 

0.8945 

0.8629 

0.8207 

0.8734 

0.8313 

0.7733 

9 

0.9026 

0.8734 

0.8345 

0.8832 

0.8442 

0.7907 

10 

0.9101 

0.8832 

0.8472 

0.8922 

0.8562 

0.8068 

Table  15 

VALUES  OF  z2  FOR  USE  IN  PREDICTIONS 
n 


For  minimum 

S'A.X ,  n 

For  maximum 

- 

dmin 

0.20 

0.26 

0.34 

0.24 

0.32 

0.43 

1 

1.504 

1  .731 

2.124 

1.650 

2.014 

2.749 

2 

1  .453 

1  .650 

1 .982 

I  .580 

1.890 

2.-91 

3 

1  .408 

1.580 

1.864 

1.520 

1  .'86 

2.284 

4 

1.369 

1  .520 

1.763 

1.467 

1.697 

2.11- 

5 

1.334 

1  .467 

1.677 

1.420 

1  .621 

1.97- 

6 

1.303 

1  .-20 

1 .60- 

1  .379 

i  .555 

1.35' 

7 

1.274 

msm 

1 .540 

1.343 

1  -*74 

1  . 7  57 

8 

1  .250 

I  .-35 

1.311 

1  .--7 

1  .672 

9 

1.227 

1  .-36 

1.282 

1  .-03 

1  .599 

10 

1.207 

1.393 

1  .256 

1.36- 

1.536 

Table  16 


PREDICTED  MAXIMUM  VALUES  OF  D 


1 

Predicted  value 

i  -  0.99  T-0-1 
(deg) 

Lower  limit 

X  -  0.96  T-0*1 
(deg) 

30.6 

29.7 

29.9 

29.0 

28 

29.3 

28.4 

29 

28.7 

27.8 

30 

28.  1 

27.2 

31 

27.5 

26.7 

32 

27.0 

26.2 

33 

26.5 

25.7 

34 

26.0 

25.3 

35 

25.6 

24.8 

36 

25.2 

24.4 

37 

24.8 

24.0 

38 

24.4 

23.6 

39 

24.0 

23.3 

40 

23.6 

22.9 

Upper  limit 

=■  1.01  T-0-1 
(deg) 


31.2 

30.5 

29.9 

29.2 

28.6 
28.  I 

27.6 

27.1 

26.6 

26.1 
25.7 

25.3 

24.9 
24.5 
24.  1 


Table  17 


DELTA  PATTERNS  GIVING  N’ON-SELF-CROSSING  EARTH-TRACKS 


26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 


■ 

n 

6 : 5 

26/13/10 

26/26/16 

26/13/4 

26/26/4 

27/27/12 

27/9/5 

27/27/19 

27/27/15 

28/14/11 

28/28/8 

28/7/2 

28/28/10 

29/29/13 

29/29/18 

29/29/6 

29/29/22 

30/15/12 

30/10/6 

30/15/5 

30/6/0 

31/31/14 

31/31/9 

31/31/22 

31/31/5 

32/16/13 

32/32/20 

32/8/3 

32/32/18 

33/33/15 

33/11/7 

33/33/7 

33/33/12 

34/17/14 

34/34/10 

34/17/6 

34/34/26 

35/35/16 

35/35/22 

35/35/25 

35/7/1 

36/18/15 

36/12/8 

36/9/4 

36/36/6 

37/37/17 

37/37/11 

37/37/8 

37/37/21 

38/19/16 

38/38/24 

38/19/7 

38/38/14 

39/39/18 

39/13/9 

39/39/28 

39/39/30 

40/20/17 

40/-0/ 12 

40/10/5 

40/8/2 

43 


Table  13 


SATELLITE  LOCATIONS  IN  24-PLANE  PATTERNS  (F$  -  2) 


Longitude  of 
ascending  node 
(all  patterns) 
(deg) 

Satellite  phase  angle  at  zero  pattern  phase  angle  (deg) 

24/24/2 

18/24/2/2.0 

18/24/2/2.667 

18/24/2/3.0 

18/24/2/7.25 

0 

0 

0 

0 

0 

0 

15.0 

30.0 

30.0 

40.0 

45.0 

108.75 

30.0 

60.0 

60.0 

80.0 

90.0 

217.5 

45.0 

90.0 

60.0 

120.0 

120.0 

120.0 

120.0 

120.0 

75.0 

150.0 

150.0 

160.0 

165.0 

228.75 

90.0 

180.0 

180.0 

200.0 

210.0 

337.5 

105.0 

210.0 

120.0 

240.0 

240.0 

240.0 

240.0 

135.0 

270.0 

280.0 

285.0 

348.75 

150.0 

300.0 

300.0 

320.0 

330.0 

97.5 

165.0 

330.0 

180.0 

0 

0 

0 

0 

0 

195.0 

30.0 

30.0 

40.0 

45.0 

108.75 

210.0 

60.0 

60.0 

30.0 

90.0 

217.5 

225.0 

90.0 

240.0 

120.0 

120.0 

120.0 

120.0 

120.d 

255.0 

150.0 

150.0 

160.0 

165.0 

228.75 

270.0 

180.0 

180.0 

200.0 

210.0 

337.5 

285.0 

210.0 

300.0 

240.0 

240.0 

240.0 

240.0 

240.0 

315.0 

270.0 

270.0 

280.0 

285.0 

330.0 

300.0 

300.0 

320.0 

330.0 

97.5 

345.0 

330.0 

°Min 

12.0° 

12.0° 

18.4° 

21.3° 

21.3° 

L_ _ 

49 


Table  19 


SATELLITE  LOCATIONS  IN  24-PLANE  PATTERNS  (F  -  4) 


Longitude  of 
ascending  node 
(all  patterns) 
(deg) 

Satellite  phase  angle  at  zero  pattern  phase  angle  (deg) 

24/24/22 

18/24/4/22.0 

18/24/4/21 .541 

18/24/4/19.667 

0 

0 

0 

0 

15.0 

330.0 

330.0 

323. 1 

295.0 

30.0 

300.0 

300.0 

286.2 

230.0 

45.0 

270.0 

60.0 

240.0 

240.0 

240.0 

240.0 

75.0 

210.0 

210.0 

203.1 

175.0 

90.0 

HT  iji 

180.0 

166.2 

1  10.0 

105.0 

120.0 

120.0 

120.0 

120.0 

120.0 

135.0 

90.0 

90.0 

83.1 

55.0 

150.0 

60.0 

60.0 

46.2 

350.0 

165.0 

30.0 

180.0 

0 

0 

0 

0 

195.0 

330.0 

330.0 

323.1 

295.0 

210.0 

300.0 

286.2 

230.0 

225.0 

240.0 

240.0 

240.0 

240.0 

255.0 

210.0 

210.0 

203.  1 

175.0 

270.0 

180.0 

166.2 

1  10.0 

285.0 

150.0 

300.0 

120.0 

120.0 

120.0 

120.0 

315.0 

90.0 

90.0 

83.1 

55.0 

330.0 

60.0 

60.0 

46.2 

350.0 

345.0 

30.0 

°Min 

21.2° 

21.2° 

o 

oo 

C4 

26.8° 
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Table  20 


SATELLITE  LOCATIONS  IN  21 -PLANE  PATTERNS 


Longitude  of 
ascending  node 

Satellite  phase  angle  at  zero  pattern  phase  angle  (deg) 

(all  patterns) 
(deg) 

21/21/19 

18/21/1/19.0 

18/21/1/18.772 

21/21/9 

18/21/0/9.0 

0 

0 

0 

0 

0 

0 

17.1 

325.7 

325.7 

321  .8 

154.3 

154.3 

34.3 

291.4 

291.4 

283.6 

308.6 

308.6 

51.4 

257.1 

257.1 

245.4 

102.9 

102.9 

68.6 

222.9 

222.9 

207.2 

257.1 

257.1 

85.7 

188.6 

188.6 

169.0 

51.4 

51  .4 

102.9 

154.3 

205.7 

120.0 

120.0 

120.0 

120.0 

0 

0 

137.1 

85.7 

85.7 

81.8 

154.3 

154.3 

154.3 

51.4 

51.4 

43.6 

308.6 

308.6 

171.4 

17.1 

17.  1 

5.4 

102.9 

102.9 

188.6 

342.9 

342.9 

327.2 

257.1 

257.  1 

205.7 

308.6 

308.6 

289.0 

51.4 

51.4 

222.9 

274.3 

205.7 

240.0 

240.0 

240.0 

240.0 

0 

0 

257.  1 

205.7 

205.7 

201  .8 

154.3 

154.3 

274.3 

171.4 

171.4 

163.6 

308.6 

308.6 

291.4 

137.  1 

137.1 

125.4 

102.9 

102.9 

308.6 

102.9 

102.9 

87.2 

257.1 

257  .1 

325.7 

68.6 

68.6 

49.0 

51.4 

51.4 

342.9 

34.3 

205.7 

24.2 
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Fig  6  Semi-log  plot  of  1/pn  v  n  for  other  tesselation* 
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Fig  10  DM||^  for  pattern  13/13/5,  at  $  *  1.0,  showing  Earth-track  in 
16-hour  orbit  (L.M  «  3:2) 
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